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ABSTRACT 

In 1870, R. Clausius found the virial theorem which amounts to introduce the trace of the 
stress tensor when studying the foundations of thermodynamics, as a way to relate the absolute 
temperature of an ideal gas to the mean kinetic energy of its molecules. 

In 1901, H. Poincare introduced a duality principle in analytical mechanics in order to study 
lagrangians invariant under the action of a Lie group of transformations. In 1909, the brothers 
E. and F. Cosserat discovered another approach for studying the same problem though using quite 
different equations. In 1916, H. Weyl considered again the same problem for the conformal group of 
transformations, obtaining at the same time the Maxwell equations and an additional specific equa¬ 
tion also involving the trace of the impulsion-energy tensor. Finally, having in mind the space-time 
formulation of electromagnetism and the Maurer-Cartan equations for Lie groups, gauge theory 
has been created by C.N. Yang and R.L. Mills in 1954 as a way to introduce in physics the differ¬ 
ential geometric methods available at that time, independently of any group action, contrary to all 
the previous approaches. 

The main purpose of this paper is to revisit the mathematical foundations of thermodynamics 
and gauge theory by using new differential geometric methods coming from the formal theory of 
systems of partial differential equations and Lie pseudogroups, mostly developped by D.C Spencer 
and coworkers around 1970. In particular, we justify and extend the virial theorem, showing that 
the Clausius/Cosserat/Maxwell/Weyl equations are nothing else but the formal adjoint of the 
Spencer operator appearing in the canonical Spencer sequence for the conformal group of space- 
time and are thus totally dependent on the group action. The duality principle also appeals to the 
formal adjoint of a linear differential operator used in differential geometry and to the extension 
modules used in homological algebra. 
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1) INTRODUCTION 


There are many ways to define the concept of ” temperature ” in thermodynamics or thermo¬ 
statics. A very useful one depends on the properties of the so-called ideal gases under a pressure 
not exceeding the atmospheric pressure, summarized by the following three experimental laws: 

1) The Boyle-Mariotte law : Discovered by Boyle in England (1662), it has been rediscovered by 
Mariotte in France (1676). For a given mass of a gas at a constant temperature, say the molar 
mass M, the product of the pressure P by the volume V occupied by this gas is (approximatively) 
constant. 

2) The Gay-Lussac-Charles law : Established around 1800 after the works of Gay-Lussac, Charles 
and Dalton, it says that, under the conditions of the preceding law, the product PV does not 
depend on the gas but only on the temperature. 

3) The Avogadro-Ampere law : Stated around 1810 by Avogadro, it says that the product PV for 
a given gas at a given temperature is proportional to the number of moles of the gas or to the num¬ 
ber of molecules as a mole is made by N molecules where N = 6, 0225.10^^ is the Avogadro number. 

As a byproduct, an ideal gas is such that PV = nRT where n is the number of moles and 
k = R/N is the Boltzmann constant while T is the ideal gas scale of temperature, also called abso¬ 
lute temperature. 

The first principle of thermostatics says that the exchange of work 6W = —PdV plus the ex¬ 
change of heat SQ = CdT -\- LdV of the system with its surroundings is a total differential, that is 
there exists a function U = U{T, V) called internal energy, such that dU = 5W -\-5Q. Accordingly, 
the properties of ideal gases are complemented by another experimental law. 

4) The Joule law : Stated by Joule in 1845 who introduced on this occasion the concept of internal 
energy, it says that the internal energy U of an ideal gas only depends on the temperature, that is 
U = U{T). 

This law has been checked by means of various expansion experiments realized by Gay-Lussac 
(1806), Joule (1845) and Him (1865). The idea is to consider an adiabatic cylinder separated in the 
middle by a wall with a tap which is suddenly opened or by a glass window which is suddenly bro¬ 
ken. One part is filled with a gas at temperature T while the other part is empty. At the end of the 
experiment, which is therefore done without any exchange of heat or work with the surroundings, 
one checks that the final temperature of the expanded gas is again T. As the new volume is twice 
the initial volume, the law follows with quite a good precision (apart for helium discovered later on). 

The second principle of thermostatics says that the 1-form 5Q admits an integrating factor 
which is a function of the absolute temperature only, that is one can find a function 9 = 9{T) and 
a function S = S{T,V) called entropy such that 8Q = 9dS has the integrating factor 1/6. More 
generally, if 6Q = 9dS = 9'dS' for two arbitrary 9{V,T) and 9'{V,T), we get S' = h{S) and thus 
1/6' = {dh{S)/dS){l/6). 

In the case of an ideal gas, dU = CdT -|- (L — P)dV dC/dV — d{L — P)/dT = 0 while 
dS = {C/6)dT + {L/6)dV ^ d{C/6)/dV - d{L/6)/dT = 0. First of all, it follows from the Joule 
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law that L = P on one side and thus C = C{T) on the other side. As a byproduct, C/9 only 
depends on T and P/9 = RT/9V must only depend on V, that is T/9 = c = cst or T = c9, a 
result showing that the ideal gas scale of temperature T can be used in place of 9 by choosing 
c = 1 or, equivalently, that the absolute temperature is only defined up to a scaling factor. It also 
follows that we may choose U = U{S, V) with dU = TdS — PdV and that the so-called free energy 
F = U — TS introduced by Helmholtz is such that dP = —PdV — SdT, a result leading therefore 
to a function F = F{T, V) allowing to define S = —dF/dT and thus U = F — TdF/dT as a way to 
bypass the principles by means of a mechanical approach to thermodynamics along the helmholtz 
analogy that we now recall. Indeed, in the lagrangian approach to analytical mechanics that we 
shall see thereafter, if one has functions q{t) of time, for example positions x{t), y{t),z{t) of points 
in cartesian space, and a lagrangian L{t,q,q) where q = dq/dt, the variational calculus applied to 
J L{t,q,q)dt may produce the Euler-Lagrange (EL) equations ~ = 0- Introducing the 

hamiltonian FI = q^ — L, we get ^ ^ and thus the conservation of energy Ff = cst whenever 

^ = 0. Accordingly, if one could find a function q{f) such that T = q, then one could recover the 
previous formulas on the condition to choose L = —F (See [22],[25] for more details). 

The following three examples are among the best ones we have been able to find in order to 
understand why exhibiting an integrating factor may not at all be as simple as what is claimed in 
most textbooks. 

EXAMPLE 1.1: {Ideal Gas) With volume V, absolute temperature T, pressure P, entropy S and 
internal energy U for one mole of a perfect gas, we obtain dU = dW + 5Q with SW = —PdV and 
6Q = CdT + LdV where C = Cy is the heat capacity at constant volume and PV = RT for one 
mole. Replacing and writing that dU and dS = {1/T)SQ are closed 1-forms, we obtain successively 
dC/dV-d{L-P)/dT = 0 and d{C/T)/dV-d{L/T)/dT = 0, that is to say A = P andC = C{T). 
We get therefore SQ = CdT -|- PdV and thus dU = CdT. However, when C is a constant as in 
the case of an ideal gas, looking for a general integrating factor of the form A{V, T)/T, the 1-form 
{CA/T)dT + {RA/V)dV must be closed and thus {C/T)dA/dV — {R/V)dA/dT = 0, a result 
leading to A = A{VT°') where a = C/R with R = Cp — Cy = (7 — ^)C according to the Mayer’s 
relation. Of course, we find the well known integrating factor 1/T leading io S = Rlog{VT°‘) and 
F = CT{1 — log{T)) — RTlog{V), but we could also use the other integrating factor VT°‘~^ leading 
to S' = RVT°‘ and get S' = R exp{S/R). If we look for an integrating factor depending only on 
T, we can only have c/T whith an arbitrary non-zero constant c used in order to fix the absolute 
temperature up to a change of scale. 

EXAMPLE 1.2:{Black Body) Using the same notations, we have now U = aVT'^ and P = 
iaT"* SQ = dU — SW = dU -\- PdV = AaVT^dT + ^aT'^dV. Looking again for any inte¬ 
grating factor of the same form A{V,T)/T as before, we should obtain AVdA/dV — TdA/dT and 
thus A = A(VT^). Of course, 1/T is the standard integrating factor leading to S' = |q;UT^ 
and F = —^aVT'^. However, VT^/T = UT^ is also an integrating factor with VT^SQ = 
AaV^T^dT + ^aVT^dV = di^aV^T'^) = dS' and we get S' = ^^S^. 

EXAMPLE 1.3: (Compare to [13], p 117) Two different ideal gases, one mole each, with re¬ 
spective heat capacities C" at constant volume V' and C” at constant volume V” such that 
C" ^ C” are contained in a cylinder and separated by an adiabatic movable piston. We shall 
prove that there cannot be any integrating factor for the exchange of heat SQ = SQ' -I- SQ” 
of this system. Using the first law of thermodynamics as in the previous examples, we have 
SQ' = C'dT' + P'dV, SQ” = C”dT” +P”dV”. However, for a reversible transformation, the piston 
must be in mechanical equilibrium and thus P = P' = P”. Now, we have PV = RT', PV” = RT” 
and we obtain therefore PdV = RdT' - {RT'/P)dP, PdV” = RdT” - {RT”/P)dP for the sys¬ 
tem now described by the only three state variables T',T”,P. Accordingly, we get the 1-form 
a = SQ = {C + R)dT' + {C” +R)dT” — {R/P){T' + T”)dP. Taking the exterior derivative, we get 
da = {R/P)dP A {dT' + dT”) and thus aAda = {R/P) {C” — C')dT' A dT” AdP ^ 0. Accordingly, 
integrating factors do not exist in general for systems which are not in thermal equilibrium. 

It remains to relate this macroscopic aspect of thermodynamics that we have presented with its 
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microscopic aspect, in particular with the kinetic theory of gases. For this, assuming the molecular 
chaos, namely that the gases are made by a juxtaposition of individual molecules of mass m with 
M = Nm, we assume that, at a given time, the directions of the speeds have a random distribution 
in space, that the size of the molecules is small compared to their respective average distance and 
that the average speed is bigger when the temperature is higher. We also assume that there is no 
interactions apart very negligible attractive forces compared to the repulsive shock forces existing 
on vey small distances. As a byproduct, the pressure is produced by the only forces acting on the 
wall of a containing volume V limited by a surface S with outside normal vector n which are made 
by the molecules hitting the surface. We explain the way followed by Clausius. 


If O is a fixed point inside V, for example the origin of a cartesian frame, and M an arbitrary 
point (for a few lines), we have: 


d{(mf 

dt 


= 2 OM . 


dOM 

dt 


2 OM .V 


d^jOMf 

di^ 


.(dOM 

^ dt ^ 


+ 2 0M . 


d^ OM 
dt"^ 


Multiplying by we recognize, in the right member, the kinetic energy of a molecule and the 

force F = acting on this molecule at time t. Summing on all the molecules contained in 

V while taking into account the fact that the sum '^m{OM)‘^ is constant when the statistical 
equilibrium is achieved, we obtain the formula: 


= ^M .F 

where the term on the right side is called virial of the gas. In the case of an ideal gas, the 
forces are annihilated two by two apart from the ones existing on S. However, the force produced 
by the pressure on a small part dS of S is known to be dF = —PftdS. Taking into account that 
P is constant inside V and on S, the total kinetic energy contained in V is thus equal to the half of: 


/ P cm .ndS = [ div{P Cm)dV = P [ div{C)M)dV = 3PV 
J S J V J V 

after using the Sokes formula because dix^ + d 2 x‘^ + dsx^ = 3. Introducing the mean quadratic 
speed u such that = Nu'^ for a mole of gas with N molecules and mass M = Nm, we obtain 

therefore PV = = ^Nmu^ = ^Mu^ and recover the experimental law found by Boyle 

and Mariotte. As a byproduct, we find = ^kT for the mean translational kinetic energy of 

a molecule. 


In order to start establishing a link between the virial theorem that we have exhibited and 
group theory, let us recall that the stress equation of continuum mechanics is drcr'"’’ = /* when the 
ambient space is with cartesian coordinates and that the stress in a liquid or a gas is the 3x3 
diagonal matrix with diagonal terms equal to —P. Using the only infinitesimal generator 6 = x’‘di 
of the dilatation group while raising or lowering the indices by means of the euclideam metric of 
R^, we obtain (Compare to (29)+(30) in [9]): 

Pdral = driPal) -a;= Pft 

as a way to exhibit the trace of the stress tensor a but, of course, it remains to justify this purely 
technical computation by means of group theoretical arguments. 

We conclude this paragraph with a few comments on the so-called axiomatic thermodynamics 
initiated by P. Duhem (1861-1916) around 1892-1894 in ([8]) and then by C. Carathodory in 1909 
([5]) (See the pedagogical review made by M. Born in 1921 [4]). 

A first comment concerns the use of differential forms (See a forthcoming paragraph for definitions), 
introduced by E. Cartan in 1899 but only used in physics and particularly in thermodynamics after 
decades. ICa = 5Q and /3 = 5W are respectively the exchange of heat and work of the system with 
its surroundings, one must never forget that any finite heat Q and work W obtained by integration 
is counted positively if it is provided to the system (One of the best references we know is by far 
[13]). In this framework, the first principle amounts io a + fi = dU where U is the internal energy 
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or, equivalently, d{a + /3) = 0. As for the second principle amounting to the existence of an ” 
integrating factor ” for a, that is the possibility to write 6Q = TdS, it is well known that it is 
equivalent to the existence of a 1-form ip = ydT such that da = p A a when n > 2 or simply 
to the condition a A da = 0 when n > 3 ([20], Th 6.4.6, p 245). Equivalently, we may use in 
both cases the Frobenius theorem saying that, for any couple of vector fields f,ri & T such that 
i{f,)a = 0,i{r])a = 0 where i() is the interior product of a vector with a form, then ii[f.,r]])a = 0 
because ii(f.)di(r]) — i{r])d(i{ff))a = i{f)i(r])da + ii[f,,r]])a,\/£,,r] G T,Va G T* from the definition 
of the exterior derivative on 1-forms. However, what is surprisingly not known at all is the link 
existing between such conditions and group theory. We start with the following key definition (See 
Section 2B and [23] for more details): 

DEFINITION 1.4: A Lie pseudogroup of transformations is a group of transformations solutions 
of a (linear or even non-linear) system of ordinary or partial differential equations called system of 
finite Lie equations. 

EXAMPLE 1.5: When n = 1 and we consider transformations y = f{x) of the real line, the 
affine group of transformations is defined by the linear system y^x = 0 with jet notations saying 
that any transformation is such that dxxf{x) = 0 while the projective group of transformations 
is defined by the nonlinear system — |(^^)^ = 0 with a similar comment. In both cases 
we have indeed a Lie group of transformations depending on a finite number of constant pa¬ 
rameters, namely y = ax + b in the first case and y = in the second case. Accordingly, 
the respective geometric object the invariance of which is characterizing the corresponding Lie 
pseudogroup are surely not made by tensors because the defining finite Lie equations are not 
first order. However, in the case of transformations of the plane (a;^,a;^) —>■ {y^,y‘^) satisfying 
y^Vi ~ y^Di = ^ J/il/i ~ 2 / 22/1 = 1 ) 110 explicit integration can be obtained 

in order to provide general solutions but another way is to say that the corresponding Lie pseu¬ 
dogroup preserves the 1-form a = x'^dx^ — x^dx^ and thus the 2-form da = dx^ A dx'^ as we have 
indeed y'^dy^ —y^dy^ = x"^dx^ — x^dx"^ and thus also dy^ Ady"^ = dx^ Adx^. The Lie pseudogroup is 
thus preserving the geometric object to = (a, da) made by a 1-form and a 2-form. More generally, 
we may consider the Lie pseudogroup preserving the geometric object uj = (a, j3) where a is a 
1 -form and /3 is a 2 -form. As da is also preserved, if we want that the system behaves at least like 
the preceding one, that is cannot have any zero order equation, we must have da = c/3 for some 
arbitrary constant c. The two pseudogroups defined by w —>■ c and oj ^ c can be exchanged by a 
change of variables bringing uj to u) if and only if c = c. This situation is the simplest example of 
the celebrated formal equivalence problem ([ 20 ],[ 21 ]). 

EXAMPLE 1.6: As a more general situation of a Lie pseudogroup of transformations of space 
with n = 3 also involving differential forms, let us consider the 1-forms a = x^dx^ and /3 = 
dx^ + x^dx^ with a A da = t),a + j3 = d(pcf + x^x^). The Lie pseudogroup preserving a and 
/3 also preserves j = da = —d(i with a A 7 = 0. It is easily seen to be made by the following 
transformations: 

2 /^ = y^ = x^ + (x^ - )x^ + a, y^ = x^/fix'^) 

where f{x^) is an arbitrary invertible function of x^ only and we have set f'{x^) = dif{x^) while 
a is an arbitrary constant, because we obtain at once -|- y^y^ = x'^ + x^x^ + a for an arbitrary 
constant a. An elementary but quite tedious computation similar to the previous one or to the ones 
that can be found in ([29], [30]) shows that solving the formal equivalence problem for w = (a,/3) 
depends on the following structure equations: 

da = Cia A P + c'l'f, djd = C 2 a A j3 + dj = c^j3 A 7 

because q;A 7 = 0, a A ft ^ 0, 7 = dx^ A dx^ ^0, /3 A 7 = dx^ A dx"^ A dx^ ^ 0. Closing this 

exterior system by taking again the exterior derivative, we get: 

Ci(ci-C3)=0, c'iC 2 - C 2 C 3 = 0 
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In the present situation, we have ci = 0, c'l = 1, C 2 = 0, C 2 = — 1, C 3 = 0. Eliminating 7 , we get the 
only conditions: 


d{a + /3) = 0, a A da = 0 

that is exactly the conditions to be found in thermodynamics through a forthcoming example. 
We invite the reader to choose: 


a = dx^, (3 = dx^ — x^dx^, 7 = dx^ A dx^ 

in order to obtain ci = = 0 , C 2 = 0 , C 2 = 1 ,C 3 = 0 with therefore no possibility to solve the 

equivalence problem u! ^ Co. 


Setting finally a = SQ = TdS,/3 = 5W with a + /3 = dU, the Helmholtz postulate, first stated 
in ([ 8 ]), assumes that it is always possible to choose the n state variables, called normal variables, 
in such a way that dT does not appear in 6W. This is a crucial assumption indeed because, intro¬ 
ducing the free energy F = U — TS, we get dF = 6W — TdS ^ S = — ^. We recall and improve 
the following result already provided in 1983 ([21], p 712-715) but never acknowledged up to now. 


THEOREM 1.7: Helmholtz postulate is a theorem whenever a A /3 ^ 0. 


Proof: Let us prove first that, setting a = T{x)dS with S' = it is always possible to choose 
the state variables in such a way that da;^ does not appear in 5W. 

Starting with n = 2, we get aA/3 = aA (a+ fi) = T{x)dx^ A dU ^ 0, implying that U must 
not involve only x^ and we may introduce the new variables = x^,y^ = U{x) in such a way 
that (3 = dy^ — T{y)dy^. Let now v(y) be a non-constant orbital integral of the ordinary dif¬ 
ferential equation dy'^/dy^ = T{y^,y‘^) satisfying therefore ^ -I- T{y)-^ = 0. It follows that 
dv = ■§^dy^ + ■§^dy‘^ = with ^ ^ 0 because otherwise we should also have ^ = 0 

and V should be constant. Using the new variables = y^,z‘^ = v(jj^,y'^), we have a Jacobian 

^ ^ 0 and obtain at once a = T{z)dz^,f3 = b{z)dz^. Finally, setting a = —SdT, we 
have da = dT A dS = da but now f3 + a = (3 + a — d{TS) = d{U — TS) = dF, that is we may 
exchange a, j3, U with a, j3, F and repeat the same procedure with T in place of S and F in place 
of U, obtaining therefore the desired result. 

Similarly, when n > 3, we can choose the new variables y^ = x^,y^ = U,y^ = x^, ■■■,y" = x" and 
obtain a = T(y)dy^,/3 = dy"^ — T(y)dy^. Considering now y^ ,..., j/" like parameters, we may use 
the same argument as above and substitute ^ = —T{y)-^ in order to get: 


. , 1 

dv = -^dy 


dv 

w 


, 2 7 3 

dy + -^dy^ 


„ dv ^ dv 


dy^ 


dv 

dy^ 


dy^ 


Choosing z^ = y^,z'^ = v{y),z^ = y^,...,z^ = y”, we obtain a = T{z)dz^,(3 = b 2 {z)dz‘^ + ... -I- 
bn{z)dz'^. The final exchange may be done as before. 

Q.E.D. 


EXAMPLE 1.8: In the case of an ideal gas with n = 2, we may choose y^ = S,y‘^ = U = CT and 
we have a = Tdy^, ft = dy"^ — Tdy^. Meanwhile, we have also TdS = CdT + PdV = dU + PdV => 
dV = ^dS - j;dU. It follows that = p, ■§^ = -p and thus = 0. Accordingly, 

V(S, U) can be chosen to be the desired orbital integral, a result highly not evident at first sight 
but explaining the notations. 

EXAMPLE 1.9: With n = 3 and local coordinates x = for the state variables, let 

us consider an abstract system with 5Q = a = x^dx^,5W = (3 = dx"^ + x^dx^. We have indeed 
da = dx^ A dx^ = ^dx^ Aa=^aAda = 0 and a + (3 = d{x‘^ + x^x^). We may therefore set 
S = x^,T = x^,U = x'^ + x^x^ and the existence of the integrating factor is compatible with the 
change of scale allowing to define T. However, we should get F = U—TS = x^ ^ dF = dx'^ and we 
should be tempted to conclude with a contradiction as we should get S = —^ = —dsF = 0 ^ x^. 
However, things are much more subtle when dealing with normal variables as it has been largely 
emphasized by Duhem in ([ 8 ]) but totally absent from the survey reference ([3]). Indeed, we 
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have now dF = 5W — SdT 5W = dF + SdT = dx^ + x^dx^ in a coherent way with the 
definition of f3. Accordingly, the correct way is thus to say that the formula S' = —^ is no 
longer true because SW now contains dT or, equivalently, that the state variables x are not nor¬ 
mal. However, exchanging U and F, it follows from our proof of the Helmholtz postulate that 
it is always possible to obtain normal state variables y = For this, we just need 

to set y^ = x^,y‘^ = + x^x^,y^ = x^ ^ x^ = y^^x^ = and obtain 

6Q = a = y^dy^,SW = (3 = dy^ — y^dy^ where (3 does not contain dy^ = dx^ = dT any longer. 
Meanwhile, dF = 5W — SdT = d{y‘^ — y^y^) = dx^ as before but now F = y"^ — y^y^ is such that 
S = x^ = y^ = = “fr as we wished. 

2) MATHEMATICAL TOOLS 

A) LIE GROUPS 

The word ''group'' has been introduced for the first time in 1830 by E. Galois and this concept 
slowly passed from algebra (groups of permutations) to geometry (groups of transformations). It 
is only in 1880 that S. Lie studied the groups of transformations depending on a finite number 
of parameters and now called Lie groups of transformations. We now describe in a modern lan¬ 
guage the procedure followed by Poincare in [19], both with the corresponding dual variational 
framework. We invite the reader to look at ([25], [26], [30], [31]) in order to discover its link with 
homological algebra and the extension functor. 

Let A be a manifold with local coordinates x = (x^,...,x") and G be a Lie group, that is 
another manifold with local coordinates a = (a^,...,a^) called parameters, with a composition 
G X G ^ G : {a,b) ^ ab, an inverse G ^ G : a ^ a~^ and an identity e € G satisfying: 

{ab)c = a{bc) = abc, aa~^ = a~^a = e, ae = ea = a, Va, b,c G G 

Then G is said to act on X if there is a map A x G —>■ A : (x,a) ^ y = ax = f{x,a) such 
that (ab)x = a{bx) = abx,'ia,b G G,Vx G A and, for simplifying the notations, we shall use 
global notations even if only local actions are existing. The action is said to be effective if 
ax = x,Vx G A a = e. A subset S' C A is said to be invariant under the action of G if 
aS d S,'ia G G and the orbit of x G A is the invariant subset Gx = {ax ] a G G} C A. If G acts 
on two manifolds A and Y, a map / : A —>■ T is said to be equivariant if /(ax) = a/(x),Vx G 
A, Va G G. For reasons that will become clear later on, it is often convenient to introduce the 
graph A x G —>■ A x A : (x, a) —>■ (x, y = ax) of the action. In the product A x A, the first factor 
is called the source while the second factor is called the target. 

We denote as usual by T = T{X) the tangent bundle of A, by T* = T*{X) the cotangent 
bundle, by A''T* the bundle of r-forms and by SqT* the bundle of q-symmetric tensors. Moreover, 
if f,ri G T are two vector fields on A, we may define their bracket [^,77] G T by the local formula 
(K.’li)* = Cdrff-’n'^dsC leading to the Jacobi identity [^, [r7,C]] + [’7, [C,C]] + [C, K,v]] = 0 ,V/,? 7 ,C G 
T allowing to define a Lie algebra. We have also the useful formula [T{f){f), T{f){r])] = T(/)([5, y]) 
where T(/) : T(A) — >• T{Y) is the tangent mapping of a map f : X ^Y. IffGT and / G G“(A), 
we set ^./ = ffdif and, if w G A''T*, we denote by i{f)uj G a'~~^T* the interior product of to by /. 
Finally, when I = {ii < ... < A} is a multi-index, we may set dx^ = dx'^^ A ... A dx*’' and introduce 
the exterior derivative d : A''T* -A a’'+^T* : uj = ojjdx^ -A duj = diOJjdx'' A dx^ with df = do d = 0 
because dijUJidx'' A dx^ A dx^ = 0 , in the Poincare sequence: 

A°T* A^T* A^T* ... A”T* —^ 0 

In order to fix the notations, we quote without any proof a few results that will be of constant 
use in the sequel (See [23] for more details). 

According to the first fundamental theorem of Lie, the orbits x = /(xo,a) satisfy the system 
of PD equations dx" jdaP = dp(x)ujl^(a) with det{uj) 0. The vector fields 9p = dp(x)9i are called 
infinitesimal generators of the action and are linearly independent over the constants when the 
action is effective. In a rough symbolic way, we have x = axo dx = daxo = daa~^x and 
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daa ^ = io = (w’’ = uj^{a)da'^) is thus a family of right invariant 1-forms on G with value in 
G = Te{G) the tangent space to G at the identity e G G, called Maurer-Cartan (MC) forms. 


Then, according to the second fundamental theorem of Lie, if 0i,...,0p are the infinitesi¬ 
mal generators of the effective action of a lie group G on X, then [0p,0^] = Cp^0T where the 
c = (Cpg. = —cfp) are the structure constants of a Lie algebra of vector fields which can be iden- 
tihed with G by using the action as we already did. Equivalently, introducing the non-degenerate 
inverse matrix a = uj~^ of right invariant vector fields on G, we obtain from crossed-derivatives the 
compatibility conditions (CC) for the previous system of partial differential (PD) equations called 
Maurer-Cartan (MC) equations, namely: 


dojf duj, 

(care to the sign used) or equivalently [ap, 




= 0 




duj^ 


2 


A = 0 




Finally, using again crossed-derivatives, we obtain the corresponding integrability conditions 
(IC) on the structure constants: 

c'^ -I- c'^ = 0 c^ c^ -\-c^ c^ =0 

^pcr ' ^crp ^5 ^pp^<7T ' '-pcT^Tp ' ^pr'-pa ^ 

also called Jacobi conditions. The Cauchy-Kowaleski theorem finally asserts that one can construct 
an analytic group G such that G = Te{G) by recovering the MC forms from the MC equations, a 
result amounting to the third fundamental theorem of Lie. 


EXAMPLE 2.A.1: Considering the affine group of transformations of the real line y = a^x-\-a^, 
the orbits are defined by a: = a^x^) of, a definition leading to dx = da^Xf^ daf and thus 
dx = {{l/a^)da^)x {daf — {a?/a})da}). We obtain therefore 0i = xdx,02 = ^ [^ 1 ,^* 2 ] = —^2 

and = {1 /a^)da^, = daf — [of /a})da} duj^ = 0, LuP' = Q ^ [ai, 02 ] = — 0:2 with 

ai = a^di of 82 , 02 = 82 . 


EXAMPLE 2.A.2: If X = a{f)xQ b{f) with a{f) a time depending orthogonal matrix (rota¬ 
tion) and b(t) a time depending vector (translation) describes the movement of a rigid body in 
K.^, then the projection of the absolute speed v = d(t)xo + b(t) in an orthogonal frame fixed in 
the body is the so-called relative speed a~^v = a~^dxo + a~^b and the kinetic energy/Lagrangian 
is a quadratic function of the 1-forms A = (a~^d, a~^b). Meanwhile, taking into account the 
preceding example, the Eulerian speed v = v(x,t) = da~^x -\- b — da~^b only depends on the 1- 
forms B = (da~^, b — da~^b). We notice that a~^d and da~^ are both 3x3 skewsymmetric time 
depending matrices that may be quite different. 


The above particular case, well known by anybody studying the analytical mechanics of rigid 
bodies, can be generalized as follows. If A is a manifold and G is a lie group (not acting necessarily 
on X now), let us consider maps a : A —>■ G : (x) —>■ (a(x)) or equivalently sections of the trivial 
(principal) bundle A x G over A, namely maps X —>■ X x G : (x) —>■ (x,a(x)). If x -|- dx is a 
point of A ’’close ” to x, then T(a) will provide a point a -I- da = a -I- fydx ’’close ” to a on 
G. We may bring a back to e on G by acting on a with a~^, either on the left or on the right, 
getting therefore a 1-form a~^da = A or daa~^ = B with value in G. As aa~^ = e we also get 
a~^da = —(da~^)a = —dbb~^ if we set b = a~^ as a way to link A with B. When there is an action 
y = ax, we have x = a~^y = by and thus dy = dax = daa~^y, a result leading to the equivalent 
formulas: 


a~^da = A = (AJ(x)dx’' = —uj'f(b(x)) 8 ib^(x)dx'') 
daa~^ = B = (Bf (x)dx^ = ujf(a(x)) 8 ia'^ (x)dx^) 

Introducing the induced bracket [A,A](f,ri) = [A(^), A(? 7 )] G G T, we may define the cur¬ 

vature 2-form dA — [A, A] = F G A^T* 0 by the local formula (care again to the sign): 


8 ^A] (x) - 8 jAl(x) - cl^Al (x)AJ(x) = Ffj(x) 



This definition can also be adapted to B by using dB + [B, B\ and we obtain: 

THEOREM 2.A.3: There is a nonlinear gauge sequence: 

X xG —^ T*^g ^ ® g 

a —^ a~^da = A —^ dA—[A,A]=F 

Choosing a ’’close” to e, that is a{x) = e + tX{x) + ... and linearizing as usual, we obtain the 
linear operator d : a°T* 0 —>■ A^T* 0 g : (A”(a;)) —>■ {diX^(x)) leading to: 

COROLLARY 2.A.4: There is a linear gauge sequence: 

A°T* ®g ^ A^T* ®g^ A^T* ®g ^ ...^ a”t* ®g —^ o 

which is the tensor product by g of the Poincare sequence: 

It remains to introduce the previous results into a variational framework. The procedure has 
been found in 1901 by H. Poincare who introduced a duality principle in analytical mechanics in 
order to study lagrangians invariant under the action of a Lie group of transformations ([19]). This 
method has been used later on by G. Birkhoff in 1954 ([2]) and V. Arnold in 1966 ([1]), each one 
omitting to quote the previous results. 

For this, we may consider a lagrangian on T* ® g, that is an action W = f w{A)dx where dx = 
dx^ A ... A dx"' and vary it. With A = a~^da = —dbb~^ we may introduce A = = —Sbb~^ G 

g = A°T*(g)C/ with local coordinates A”(a;) = —u}"{b{x))Sb'^{x) and we obtain 6A = dX—[A, X] that 
is SAJ = diX" — c"^A^X‘^ in local coordinates. Then, setting dw/dA = A = {A\.) G a”“^T* 0 g, 
we get: 


SW = j A6Adx = j AidX - [A, X])dx 

and, after integration by part, the Euler-Lagrange (EL) relative equations ([22],[23]): 

+ c;,APA^ = o 

Such a linear operator for A has non constant coefficients linearly depending on A and the struc¬ 
ture constants. Setting Saa~^ = /i G t/, we get A = a~^{Saa~^)a = Ad{a)^ while, setting 
a ^ a' = ab, we get the gauge transformation A A' = (ab)~^d(ab) = b~^a~^(dab + adb) = 

Ad(b)A + b~^db,\/b G G. Setting b = e + tA + ... with t <C 1, then 5A becomes an infinitesimal 

gauge transformation. However, setting now a ^ a' = ca, we get A' = a~^c~^{dca + cda) = 

a~^{c~^dc)a + A and thus 5A = Ad{a)dfj, when c = aba~^ = e + t/x + ... with t <C 1 (See [23], 

p 180, 424 for more details and compntations using local coordinates). We may also notice that 
aa~^ = e + a(l(a“^) = 0 ^ 6{a~^) = —a~^6aa~^ and thus: 

5A = —a“^5aa“^(ia + a“^(i(((laa“^)a) 

= —a~^Saa~^da + a“^(i((laa“^)a + a~^6aa~^da 
= Ad{a)dpL 

Therefore, introducing by duality B such that Bp, = AX, we get the divergence-like absolute equa¬ 
tions diB]j = 0. When n = 1, we recognize at once the Birkhoff-Arnold dynamics of a rigid body, 
with time t as independent variable, or the Kirchhoff-Love theory of a thin elastic beam, with 
curvilinear abcissa s along the beam as independent variable. 

REMARK 2.A.5: As the passage from A to B, that is from left invariance to right invariance is 
not easy to achieve in actual practice, we indicate a way to simplify the use of the adjoint mapping 
(Compare to [23], Proposition 10, p 180). Indeed, working formally, from SA = dX — cAA, we 
may define on G a square matrix acting on g and define p hy X = Mp. Substituting, we obtain 
6A = d{Mp) — cAMp = Mdp -\- {dM — cAM)p and thus 6A = Mdp AA dM — cAM = 0, that is to 
say M{b) must be a solution of the linear system of PD equations -I- c"p^ujP{b)M^ = 0. It just 
remains to prove that this system is involutive by computing the crossed derivatives. An easy but 
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tedious computation provides: 


db^ 






)-(< 


'per ^a(3 


■ CpaC^^ + c"ppCaa)<^>s)M'^ = 0 


Using both the MC equations and the Jacobi conditions achieves the proof of this technical but 
quite useful result. 


We may therefore ask: 

PROBLEM: HOW IS IT POSSIBLE AND WHY IS IT EVEN NECESSARY TO INTRODUCE 
DIEEERENT EQUATIONS WITHIN THE SAME GROUP BACKGROUND. 

B) LIE PSEUDOGROUPS 

We start recalling a few notations and definitions about fibered manifolds and their jet bundles 
(See [20] and [22] for more details). In particular, if f —>• A : (x, y) —>■ (x) is a fibered manifold with 
changes of local coordinates having the form x = (f{x), y = 'tp{x, y), we shall denote by Jq{£) —>■ X : 
{x,yq) —>■ (x) the q-jet bundle of S with local coordinates for i = 1, ...,n, k = 1, 0 <| 

y.\< q and y^ = y’^. We may consider sections fq : (x) (x, f^(x), f,^(x), f^^/x),...) = (x,fq(x)) 

transforming like the sections jq(f) : (x) —>■ (x, f'^(x), dif'^(x), dij f^(x),...) = (x,jq(f)(x)) where 
both fq and jq{f) are over the section / : (x) —>■ {x,y^ = f^{x)) = {x,f{x)) of S transforming 
like f{(p{x)) = 'tp{x,f{x)). If T{S) has local coordinates {x,y]u,v), we shall denote by V(£) the 
vertical bundle of £, namely the sub-vector bundle of T{£) with local coordintates (x, y; 0, v). The 
(nonlinear) Spencer operator just allows to distinguish a section fq from a section jq{f) by in¬ 
troducing a kind of ” difference ” through the operator D : Jq+i{£) —>• T* 0 V{Jq{£)) : fq+i 
jiifq) - fq +1 with local components ((9*/'=(x) - /f (x), 9i/j=(x) - /^(x),...) and more generally 
(Dfq^i)^ -(x) = 9i/^(x) — /^_|_i.(x). If m = n and £ = X x X with source projection, we denote 
by Hq = \iq{X,X) C Jq{X X X) the open sub-bundle locally defined by det{yf) ^ 0 and we shall 
set A = det{dif^{x)). Also, if £ and A are two fibered manifolds over X with local coordinates 
(x,j/) and (x,z) respectively, we shall denote by £xxX their fibered product over X with local 
coordinates {x,y,z). Finally, if E is a vector bundle over X with transition rules having the form 
X = (p{x), y = A(x)y, we shall denote by E* the vector bundle obtained from E by inverting the 
transition matrices, exactly like T* is obtained from T. 

In 1890, Lie discovered that Lie groups of transformations were examples of Lie pseudogroups 
of transformations along the following definition which expands the preliminary Definition 1.4: 

DEFINITION 2.B.1: A Lie pseudogroup of transformations P C aut{X) is a group of trans¬ 
formations solutions of a system of OD or PD equations such that, if y = /(x) and z = g{y) are 
two solutions, called finite transformations, that can be composed, then z = g o f(x) = h(x) and 
X = f~^{y) = g(y) are also solutions while y = x is the identity solution denoted by id = idx 
and we shall set ida = ja(id). In all the sequel we shall suppose that P is transitive that is 
Vx,j/G A,3/er,y = /(x) 

We notice that an action y = /(x, a) provides a Lie pseudogroup by eliminating the p param¬ 
eters a among the equations yq = jq(f)(x,a) obtained by successive differentiations with respect 
to X only when q is large enough. The system TZq C H, of OD or PD equations thus obtained 
may be quite nonlinear and of high order. The concept of parameters is not existing in this new 
framework and thus no one of the methods already presented may be used any longer. Setting 
fix) = /(x, a(x)) and /q(x) = jq(/)(x, a(x)), we obtain a(x) = a = cst fq = jqif) because 
Dfq^i = jiifq) — fq +1 = (9/q(x,a(x))9i(x)(ix* as a 1-form and the matrix involved has 
rank p in the following commutative diagram: 

TZq ^0 

Jqif) nt fq 
X 


X xG 

a = cst tit O'ix) 
X 
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More generally, looking now for transformations ’’close” to the identity, that is setting y = 
X + t^(x) + ... when t <C 1 is a small constant parameter and passing to the limit t — >■ 0, we 
may linearize any (nonlinear) system of finite Lie equations in order to obtain a (linear) system 
of infinitesimal Lie equations Rq C Jq(T) for vector fields. Such a system has the property that, 
if T] are two solutions, then [^, 77 ] is also a solution. Accordingly, the set 6 C T of its solutions 
satisfies [0,0] C 0 and can therefore be considered as the Lie algebra of L. 

Looking at the way a vector field and its derivatives are transformed under any / G aut{X) 
while replacing jq(f) by fq, we obtain: 

= frix)C{x) ^ vtifix))ffi{x) = fl:{x)f^{x) + fl:i{x)Cix) 

and so on, a result leading to: 


LEMMA 2.B.2: Jq{T) is associated with Ilg+i = ng+i(A, X) that is we can obtain a new section 
rjq = fq+i{fq) from any section fq G Jq(T) and any section /g+i G n^+i by the formula: 

+ ... = frX + ... + <\y\<q 

where the left member belongs to V{Ilq). Similarly Rq C Jq{T) is associated with TZq+i C n,+i. 


We now need a few basic definitions on Lie groupoids and Lie algebroids that will become 
substitutes for Lie groups and Lie algebras. The first idea is to use the chain rule for derivatives 
jq(.9 o /) = jqid) o jq{f) whenever f,gG aut(X) can be composed and to replace both jq{f) and 
jq{g) respectively by fq and gq in order to obtain the new section gq o fq. This kind of ’’compo¬ 
sition” law can be written in a pointwise symbolic way by introducing another copy Z of X with 
local coordinates (z) as follows: 


Oz Ov dz dv 

Iq ■ Uq{Y,Z)XYUq{X,Y) -)> Uq{X, Z) '. {{y,z, ^)...)> ix,y, —,...) {x,z, 

We may also define jq{f)~^ = jq{f~^) and obtain similarly an ’’inversion” law. 


DEFINITION 2.B.3: A fibered submanifold TZq C Ilq is called a system of finite Lie equa¬ 
tions or a Lie groupoid of order q if we have an induced source projection Uq : TZq -G X, target 
projection j5q : TZq -G X, composition 7 ^ : TZqXxTZq — >■ TZq, inversion iq : TZq — >■ TZq and iden¬ 
tity idq : X —>■ TZq. In the sequel we shall only consider transitive Lie groupoids such that 
the map {oiq,(3q) : TZq —>■ A x A is an epimorphism. One can prove that the new system 
Pr{TZq) = TZq+r = Jr (TZq) fl Tlq+r C Jrijfq) obtained by differentiating r times all the defining 
equations of TZq is a Lie groupoid of order g + r. 

Using the algebraic bracket {jq+i{f), jq+ifr])} = jq{[£,,y]),'Jf,y G T, we may obtain by bilin¬ 
earity a differential bracket on Jq{T) extending the bracket on T: 

i^q^Vq] = {^q-el,Vq-el} + 'iiODyq+l - i{y)Dfq+l ,yfq, yq G Jq{T) 

which does not depend on the respective lifts and 77^+1 of fq and yq in Jq^iiT). One can 
prove that this bracket on sections satisfies the Jacobi identity and we set: 


DEFINITION 2.B.4: We say that a vector subbundle Rq C Jq{T) is a system of infinitesimal 
Lie equations or a Lie algebroid if [Rq,Rq] C Rq, that is to say [fq,yq] G Rq,'ifq,yq G Rq. Such a 
definition can be tested by means of computer algebra. We shall also say that Rq is transitive if we 

tt’ 

have the short exact sequence 0 —>■ —>■ A T —>• 0. In that case, a splitting of this sequence, 

namely a map Xq '■ T ^ Rq such that o Xq = or equivalently a section Xg G T* 0 Rq 
over idr G T* Z) T, is called a Rq-connection and its curvature Kq G A^T* Z Rq is defined by 
Kqif,y) = [XqiOXqiv)] “ Xq{[^,V]),'J^,V T. 

PROPOSITION 2.B.5: If [Rq,Rq\ C Rq, then [Rq+^^Rq+r] C i?g+r,Vr > 0. 
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Proof: When r = 1 , we have pi{Rq) = Rq+i = {C9+1 G Jq+i{T) \ ^ G Rq,DS^q+i & T* ® Rq} 
and we just need to use the two following formulas showing how the Spencer operator acts on the 
various brackets (See [GBl] [deform] for more details): 

iiODi^q+l^Pq+l} = {i{C)D^q+l,'qq} + {S^q^i^QDpq+l}, \/C,&T 

*(C)^[C9+l,??9+l] = [iiOD^q+l,Vq] + [^q, iiC)Dpq+l] + i{L{r]l)C)D^q+l - Dpq+l 

because the right member of the second formula is a section of Rq whenever 1,775+1 € Rq+i- 
The first formula may be used when Rq is formally integrable. 

Q.E.D. 

EXAMPLE 2 .B. 6 : With n = l,q = 3 ,X = R and evident notations, the components of [$3,773] 
at order zero, one and two are defined by the totally unusual successive formulas: 

[^,v] = ^d:,r] - 

= ^9x1]x V^x^x 

{l^2t'n2\)xx — ^xVxx ‘Hx^xx “t“ ^dxT]xx ‘H^x^xx 

([6,^3]) XXX - ‘^^x'Hxxx x^xxx + $ 9 x 77 

XXX - ??9x$ XXX 

For affine transformations, $xx = 0,77^^ = 0 => {[^2,V2\)xx = 0 and thus [-R2,^2] C i?2- 

For projective transformations, ^xxx = 0 ,r]xxx = 0 i[^3,'n3])xxx = 0 and thus [i?3,-R3] C R3. 

The next definition will generalize the definition of the classical Lie derivative: 

£($)w = ( 7 ($)d + d 7 ($))w = ^jqiexp t$)“^(w)lt=o- 

DEFINITION 2 .B. 7 : We say that a vector bundle F is associated with Rq if there exists a first 
order differential operator +($5) : + —>■ + called formal Lie derivative and such that: 

1) '9($5 + rjq) = L{fq) + L{riq) '^^qiVq G Rq- 

2 ) L{fiq)=fL{iq) yfqGRqyfGC^iX). 

3 ) [- 9 ($g),L(? 75 )] =L($ 5 )oL( 775 )-L(? 75 )oL($ 5 ) =+([$5,77,]) ^Jfq^Pq & Rq. 

4 ) +($,)(/??) = fL{fq)v + {U)v V$5 G i?5,V/ G C°°(A),Vr7 G F. 

As a byproduct, if E and F are associated with Rq, we may set on E ^ F: 

L{fq)iv'S)C) = L{fq)ri'S) C + L{^q)C V$5 G Rq,'i'q & E,\/C G F 

If 0 C r denotes the solutions of Rq, then we may set £($) = L(j9($)),V$ G 0 but no explicit 
computation can be done when 0 is infinite dimensional. 

PROPOSITION 2 .B. 8 : Jq{T) is associated with J5+i(T) if we define: 

L {(, q + i)'nq = {$5+1,775+1} +i($)D? 75 +i = [$5,775] +i( 77 )D$ 5 +l 

and thus Rq is associated with R5+1. 

Proof: It is easy to check the properties I, 2 , 4 and it only remains to prove property 3 as follows. 

[' 9 ($<?+l); - 9 (’^<?+l)]C<? “ 9/($5 +i)({775+i, $5+1} + i{r])Dfq.^.l) — Z/(?75+i)({$5+i, $5+1} + 7($)D$5+i) 

{$<7+11 {^<7+2 ; C(/+2 }} {^(7+1 7 {$<7+25 $(7+2}} 

+ {^q+l,'i{v)DCq+2} - {775+1, i($)DC5+2} 

+ 7 ($)D{? 75 + 2 , C9+2} - i( 77 )D{$ 5 + 2 , Cg+2} 

+i($)D(i( 77 )DC 5 + 2 ) - 7 (? 7 )D(i($)DC 5 + 2 ) 

“ {{$<7+2) ''7<7+2}) C<7+l} + {*($)4^''7<7+2) C<7+l} ~ {'liv)Df,q+2t C<7+l} 

+'i([^,v])DCq+i 

~ {[$ 9 + 1 ) l?<J+l]) C< 7 +l} + *([$) 17 ])^C< 7+1 

by using successively the Jacobi identity for the algebraic bracket and the last proposition. 

Q.E.D. 
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EXAMPLE 2.B.9: T and T* both with any tensor bundle are associated with Ji{T). For T we 
may define = [C,??] + = HiJiiv)}- We have = 

~V^^s + Cdrrf and the four properties of the formal Lie derivative can be checked directly. Of 
course, we find back £{^)ri = [C,? 7 ],V ^,?7 G T. We let the reader treat similarly the case of T*. 


THEOREM 2.B.10 : There is a first nonlinear Spencer sequence: 

0 ^ aut{X) n,+i(A, A) A T* (g) JfiT) A 0 Jg-i(T) 
which is locally exact if A ^ 0, with restriction: 

0 ^ r A Uq+i A T* (g) i?,j A A^T* (g) J,_i(r) 

Proof: There is a canonical inclusion Ilq+i C Ji(ng) defined by j/* j = and the composition 

ffi+i ° ji{fq) is a well defined section of Ji(nq) over the section f~^ ° fq = idq of 11^ like idq+i. 
The difference Xq = ^/<?+i = ffi+i ° Ji(/g) ~ idq+i is thus a section of T* (g) ^(11,) over idq and 
thus of T* (g) Jq(T). For <7=1, setting xo = A — id € T* 'S) T and gi = we get: 

x’} = gfd^f - = Al - St xli = g’tml - au^) 

We shall prove later on the useful formula frX'p.y + + fjl+irX^i = ^i//J ~ Z^+ii allowing to de¬ 

termine Xq inductively. 

We have D'Xqi^^v) = Dxq{^,g) — {XqiO^Xqiv)} = 0 and provide the only formulas that will be 
used later on and can be checked directly by the reader: 



- djx'^i - xlj + Xj,i - 
^jxti ~ Xu,j + X%,i — {x'",iXw,j 


- {XyXr,j - X^jXr,i) = 0 
+ xiiXrj - xljXr,i - Xjxty) = 0 


( 1 ) 

( 2 ) 


In these sequences, the kernels are taken with respect to the zero section of the vector bundles 
involved. We finally notice that the condition det{A) 0 amounts to A = det{dift 0 because 
det{fl^) 0 by assumption. One can prove by induction that the first nonlinear Spencer sequence 
is locally exact if det{A) 0, that is any section of T* (g) Jq(T) killed by D' is locally the im¬ 
age by i) of a section of n^+i, contrary to its restriction (See [23], p 215 for more details and 
compare to [14], p 162, 195). Also, introducing the vector bundle Ci = T* ® Rq/Sgqj^i, we have 
det(A) 0 3xq = TqO A with ^ = X^n,r{A~t\ = ~gifl+ii + and D induces a nonlinear 
operator Di : —>■ Ci, a result that will be generalized later on in the linear framework. The 

brothers Cosserat were speaking about the lagrangian field Xq and the eulerian field Tq defined in 
([7], §71, (70)-|-(71) O (72)-|-(73), p 190). This is a subtle confusion because the true eulerian field 
Uq = —Dfqti^ obtained by exchanging source with target, cannot be expressed from Xq by means 
of linear algebra (See [22], p 303 for more details). 

Q.E.D. 


REMARK 2.B.11: Rewriting the previous formulas with A instead of xo we get: 


d^A] - djAl - Al 




= 0 


~ djxt ~ Xl,iXr,j + X\,jXr,i ~ A'lxwj + A^jXw,i — b 


( 1 *) 

( 2 *) 


When <7 = 1 and 772 = 0, we find back exactly all the formulas presented by E. and F. Cosserat in 
[22], p 123 and [34]) (Compare to [14]). We finally notice that Xq = ~Xq is a i?q-connection if and 
only if A = 0, a result in contradiction with the use of connections in physics (Compare to [14], p 
162, 195). However, when A = 0, we have XoiO = C and thus (exercise): 


D'Xq+l = iDXq+l)i^,V) - {[Xq{0,Xq{g)] + i{0D{Xq+l{v)) - i{g)D{Xq+liO)) 
= -[XqiO,Xqiv)]-Xqi[^,V]) 

does not depend on the lift of Xq- 
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THEOREM 2.B.12: In the case of a lie group of transformations, the nonlinear Spencer sequence 
is isomorphic to the nonlinear gauge sequence when q is large enough and we have the following 
commutative diagram ([22], [23]): 


X xG 

-+ 

T*®G 

MC 

-+ 

XT* ® G 

i 




i 

TZq 

D 

-+ 

T* ®Rq 

D' 

-+ 

XT* ® Rq 


0 r - 

The action is essential in the Spencer sequence but disappears in the gauge sequence. 


Introducing now the Lie algebra G = 7e(G) and the Lie algebroid Rq C Jq{T), namely the 
linearization of TZq at the q-jet of the identity y = x, we get the commutative and exact diagram: 


0-J> X X G = Rq ->0 

A = cstntA(a:) J,(e)nte9 

X = X 

where the upper isomorphism is described by —>■ = X{x)d^9’^{x) for q large enough. 

Applying the Spencer operator, we finally obtain di^’^{x) — (x) = diX{x)dfj,9^{x) and get: 

COROLLARY 2.B.13: The linear Spencer sequence is isomorphic to the tensor product of the 
Poincare sequence by G in the following commutative diagram: 

Xt*®G — Xt*®G — Xt*®G 

'I' 

0^0^ XT* ®Rq ®Rq XT* ® Rq 

where the vertical isomorphisms are induced by the previous diagram. 

When E is a vector bundle over X and D : Jqj^i{E) —>• T*®Jq{E) is the corresponding (linear) 
Spencer operator, we denote by 5 : S'^+iT* ® E ^ T* ® SqT* 0 E the Spencer map induced (up 
to sign) by applying D to the short exact sequence 0 —>■ Sq+iT* 0 E —>• Jg+i(E) —>• Jq{E) —>• 0. 
We can extend the Spencer operator to an operator D : XT* ® Jg+i(E) —>• a’’+^T* 0 Jq{E) : 
a(8)^q+i —>■ da®^q-\-{—\Ya/\D^q+i and the corresponding Spencer map 5 : XT* ^^Sq+iT* ®E 
A^’+ip* 0 SqT* (g) E is defined by (dw)^ = dx* A . For any linear system Rq C Jq{E), we 

may define the r-prolongation pr{Rq) = Rq+r = Jr{Rq) H Jq+r{E) C Jr{Jq{E)) and the symbol 
Prigq) = gq+r = Rq+r H Sq+X* ® E both with the restrictious D : XT* 0 Rq+i —>■ a’'+^T* 0 Rq 
and d : XT* ® —>• a’'+^T* ® gq. It is finally easy to verify that = 0 ^ = 0 ([20], [25]). 

DEFINITION 2.B.14: A system Rq C Jq{E) is said to be formally integrable if all the equa¬ 

tions of order q + r that can be obtained from the system are obtained by differentiating r times 
only the equations of order q defining Rq or, equivalently, if the maps Rq+r+i Rq+r are epi- 
morphisms Vr > 0. Its symbol gq C SqT* 0 E is said to be finite type if gq+r = 0 for r large 

enough, I-acyclic if all the sequences ... -A A‘^T* 0 gq+r -A ... are exact Vr > 0,Vs = l,...,l and 

involutive if it is n-acyclic. A finite type symbol gq is involutive if and only if pq = 0. Finally, 
a system is said to be involutive if it is formally integrable and if its symbol is involutive. Such 
crucial properties can now be checked by means of computer algebra techniques based on the 
Janet/Goldschmidt/Spencer criterion saying roughly that Rq is formally integrable whenever pq is 
involutive or even 2-acyclic and : Rq+i -A Rq is an epimorphism. Otherwise, one may start 

afresh with Rq^^ = 7r^+^(Eq+i) C Rq and so on, till the criterion could be used ([12], [20], [25]). 

EXAMPLE 2.B.15: Linearizing the finite Lie equations of Example ..., we find a system 
El C Ji{T) defined by the two first order equations x'^fl — x^ff -I- = 0,x^^2 ~ ~ = 0- 

In such an example, pi is involutive (exercise) but the system is not formally integrable because, 
using crossed derivatives, one can obtain the new first order equation ^ 1 +^ 2 = 0. The combined 
first order system, namely the projection of E 2 into Ei, is involutive with the same solutions. 
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REMARK 2.B.16; The (first) linear Spencer sequence 0 —>■ E Jq+i{E) T* 0 Jq{E) 

/\^T* ® Jq-i{E) and its restriction 0 —>■ 0 Rq+i T*®Rq /\^T* ® Jq-i{E) are not very 
useful in actual practice because the operator D is not involutive and even not formally integrable. 
Indeed, from the first order equations = 0, we obtain, by using crossed derivatives, the 

new first order equations di^j — = 0. 


For any involutive system Rq C Jq{E) the Janet bundles = f\''T* ® Jq{E)/ ® Rq + 
(5(A’'“^T* 0 Sq+iT* 0 E)) and Spencer bundles Cr = A’’T* 0 Rq/S{A'^~^T* 0 Pq+i) C Cr(E) = 
A'~T* (g) Jq{E)/S{A'~~^T* (g) Sq+iT* 0 E) are related by the short exact sequences 0 ^ Cr —>■ 

Cr{E) -5^ FV —>■ 0 where the epimorphisms d’r for r = 0,1, ...,n are induced by the epimorphism 
$ = $0 : Co{E) = Jq{E) — >• Jq{E)/Rq = Fq. In the commutative diagram below where all the 
operators are induced by D, the (second) linear Spencer sequence is the kernel of the projection of 
the linear hybrid sequence onto the linear Janet sequence'. 


0 ^ 


0 0 


0 






i 


0 


0 


Co 

A 





i 


0 


E 


CoiE) 

A 



II 


i^o 


0 


E 


Fo 

A 


i 

0 





i 


Cl 

A 


Cn 

^ 0 

i 



i 


Ci{E) 

A 


CniE) 

^ 0 




i 


El 

A 


Fn 

^ 0 

i 



i 


0 



0 



In particular, if E = T and Rq C Jq{T) is a transitive involutive system of infinitesimal Lie equa¬ 
tions, the Janet bundles are associated with Rq. If moreover pq = 0, then, whenever the dimension 
of the underlying Lie group is increasing, the dimensions of the Janet bundles are decreasing while 
the dimensions of the Spencer bundles are increasing by the same amount. We obtain therefore 
the following picture: 


SPENCER 

\ 

o 

II \ 

_ II JANET 

showing why, in some virtual sense, Janet and Spencer are playing at see-saw ([12]-0'[33]). This 
picture will give the key for all the applications we shall present in the next section. 

EXAMPLE 2.B.17: When n = 3 and E = X xM., the second order system R 2 C J 2 {E) defined 
by the three PD equations 1/33 = 0,j/23 — yii = 0,7/22 = 0 is trivially formally integrable because 
it is homogeneous but is not involutive because its symbol g 2 with dim{g 2 ) = 6 — 3 = 3 is finite 
type with dim{g^) = 1 and g^+r = 0, Vr > 0. Accordingly, we have dim{R 2 ) = l-l-3-|-3 = 7 while 
dim{Rz+r) = 8 = 2”,Vr > 0 ([16], p 79). We let the reader prove as an exercise of linear algebra 
that 53 is 2-acyclic by showing the exactness of the J-sequence 0 —>■ A^T* 0 <73 -A A^T* 0 (/2 —>■ 0 
and we may consider the first prolongation i?3 C JsiE) defined by the following 12 PD equations: 


{ = 2/333 = 0 

1 2 3 

1 4 >^ = 2/233 = 0, (/>^ = 2/223 = 0, c/)"* = 2/222 = 0 

1 2 • 

1 = 2/133 = 0, = 2/123 - 2/111 = 0, <()’^ = 2/122 = 0, = 7/113 = 0, = 2/112 = 0 

1 • • 

[ = 2/33 = 0, = 7/23 - 2/11 = 0, = 7/22 = 0 

• • • 


In this particular situation, that is when g^ is already 2-acyclic though NOT involutive, it is known 
that the generating compatibility conditions (CC) are first order (See [23], p 120) and described by 
the following 21 PD equations: 
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= dzpP — da^^ = 0,..., = (i30® — -I- di^®,..., = d^cp^ — da(/)® = 0 

1 2 3 

= da<()® — di(l)2 = 0 , dtio = da0® — di^® -|- dicp^ = 0 ,..., = daf/)® — di(jf^ = 0 

1 2 • 

= d3^®® — (p^ = 0, ip^^ = d^cp^^ — (p"^, di*®® = d3(()®® — (p^ = 0 

1 2 3 

= da4>®° - = 0 , = d2())“ - -b = 0 ,= d2(p^^ - $4 = 0 

I 2 • 

-i/)^® = di^®® — (p^ = 0, 'i/’®® = di^®® — (p^ = 0, = di(()®® —4^ = 0 

I • • 


Each dot is producing one CC apart from one as we may verify the relation: 


= d22</>® - 

and check therefore the remaining 13 — 1 = 12 first order CC: 


f d® = dstp^ — da'i/'"* + dif/'® = 0,..., d® = dstp^^ — daV'® + diip^^ = 0 

I 2 3 

1 d^ = d3V’^® - daV-^^ + 4^ =0 ,..., d® = d3V’^® - daV’^® + 4^ = 0 

I 2 3 

1 d”^ = d3^®® — di4^^ + 4"^ = ■■■, d® = d3^®® — di4^^ -b ^® = 0 

I 2 3 

[ d®® = daV’^® — dif/'^® +4^ = d®^ = daV"^^ — di4^^ + 4^^ = 0 

I 2 • 


It is quite a pure chance that this system is involutive with the following 3 first order CC: 

{ dsO^^ - + did^ -6^=0, dgd^^ - dad® + did® - d^ = 0 I 1 2 3 

THE FOLLOWING ABSOLUTELY NONTRIVIAL POINT WILL BE CRUCIAL FOR UNDER¬ 
STANDING THE STRUCTURE OF THE CONFORMAL LIE EQUATIONS LATER ON. 
Indeed, with q = 3 and 34 = 0, we can define the Spencer bundles to be Cr = A’T* (g)i? 3 , construct 
in any case the Janet sequence for the trivially involutive operator J 3 and obtain the following 
contradictory diagram where dim{F^) = 2 instead of the awaited 3 (!): 






0 


0 


0 


0 






i 


i 


i 


i 


0 


0 

4 

8 

4 

24 

4 

24 

4 

8 

^ 0 





i 


i 


i 


i 


0 


I 

4 

20 

4 

45 

4 

36 

4 

10 

^ 0 



II 


i$o 


i $1 


i $2 




0 


I 

V 
—)■ 

12 

4 

21 

4 

12 

4 

2 

^ 0 





i 


i 


i 


i 






0 


0 


0 


0 



The explanation needs difficult homological algebra even on this elementary example which could 
be nevertheless treated by means of computer algebra while using quite large matrices. Indeed, 
starting from the short exact sequence 0 —>■ R 3 —>■ Jz{E) Fq —>■ 0 with fiber dimensions 
0—>- 8 —>-20—>-12—>-0 and using 3 prolongations in order to ’’reach” F 3 , we get the following jet 
sequence of vector bundles, in fact the same that should be produced by any symbolic package: 

0 —>■ i?6 Je{E) —>■ J 3 {Eo) —>■ J 2 {Ei) —>■ Ji{F 2 ) —>■ F3 —>■ 0 
with respective fiber dimensions: 

0 8 84 -)> 240 210 48 -)> dimiFs) 0 

Accordingly, if the sequence were exact, using the Euler-Poincare formula ([15], Lemma 2.2, p 206), 
we should get dimlFs) = 48 — 210 -I- 240 — 84 -|- 8 = 2, a result showing that the sequence cannot 
he exact. Knowing why it is not exact and what is the resulting cohomology needs the following 
diagram obtained by induction, where all the rows are exact but perhaps the upper one: 
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0 

0 

0 

0 


0 


i 

i 

i 

i 


i 

0 ^ 

56 

^ SqT* 

^ S3T* (S Eo 

^ S2T* (S Fi 

-A 

T* SF2 


id 

id 

id 

id 


II 

0 ^ 

T* ®gs 

T*(S S5T* 

T* (g) S2T* (g) +0 

T* ST* SFi 

-A 

T* (g) +2 


id 

id 

id 

id 


i 

0 ^ 

A^T* ® 54 

^ A^T* (S SiT* 

A^T* (g) T* (g) +0 

-)> A^T* (g) Fi 

-A 

0 


id 

id 

id 

i 



0 

A^T* (g) 53 

-)> A^T* (g) S3T* 

A^T* (g) Fo 

0 




i 

i 

i 





0 

0 

0 




As 54 

= 95 = ge 

= 0 and dim{A^T* 

(g) 53 ) = dimigs) = 1, 

a chase using the standard snake 


lemma of homological algebra ([32], p 174) proves that the upper seqence is not exact at S 2 T* (Si Fi 
with cohomology of dimension 1. Hence, the previous sequence is not exact at J 2 {Fi), that is with 
dim{im{J 3 {Fo) —>• J 2 {Fi))) = 240—84+8 = 164 while dim{ker{J 2 {Fi) —>• Ji(^ 2 ))) = 164+1 = 165 
and we have indeed 48 — 210 + 165 = 3. 

The explanation of this tricky situation is not easy to grasp by somebody not familiar with 
homological algebra. Indeed, let us apply the d-map inductively to the short exact sequence 
0 —>■ gq+r —t Sqj^rT* (S E ^ hr ^ 0 and consider the right part of the diagram thus obtained 
where the middle row is exact (See [23], p 151,152 for more details): 

0 0 

... ^ (S Sq+iT* (S E 4 A^T* (S SqT* (S E ^0 

i ; 

... ^ (Shi 4 A”r* (g) +0 

i 

0 

Cutting the diagram, we may consider the following quotient diagram: 

0 0 


0 ^ 

d(A"-iT*®g,+i) 


A”r* (g) Rq 

-)> 

Cn 

^ 0 


i 


i 


i 


0 ^ 

d(A”-iT* (g) Sq+iT* (g) E) 


A”r* (g) Jq{E) 

-)> 

Cn{E) 

^ 0 


i 


i 


i 


0 ^ 

d(A”-iT* (g)/ii) 


A”T* (g) +0 


Fn 

^ 0 


i 


i 


i 



0 


0 


0 



When gq is n— 1-acyclic but NOT n-acyclic, then hi is NOT n — 1-acyclic and a chase is showing 
that the left vertical column is not exact at the central vector bundle. Using again the snake 
lemma, there is no way to get an upper injective map in the right vertical column. In the present 
situation with n = 3 and g = 3, as gs is 2-acyclic but NOT 3-acyclic and 34 = 0, we have indeed 
dim{F^) = 3 because dim(5{A^T*(Shi)) = 10—1 = 9 in a coherent way with explicit computations. 
Accordingly, the only correct diagram allowing to deal with exact sequences on the jet level is the 
following one where all the operators involved are involutive, the n = 4 vertical sequences are short 
exact sequences and 1 — 27 + 60 — 46 + 12 = 8 — 24 + 24 — 8 = 0 ([15], Lemma 2.2, p 206): 
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0 


0 


0 


0 






i 


i 


i 




0 


e 

4 

8 

4 

24 

4 

24 

4 

8 

^ 0 





i 


i 


i 


i 


0 


1 

4 

35 

4 

84 

4 

70 

4 

20 

^ 0 



II 


i$o 


i $1 


i 4*2 




0 ^ 0 


1 

V 
—)■ 

27 

4 

60 

4 

46 

4 

12 

^ 0 





i 


i 


i 


i 






0 


0 


0 


0 



As a byproduct, ONE MUST CONSTRUCT THE JANET AND SPENCER SEQUENCES FOR 
AN INVOLUTIVE SYSTEM IN ORDER TO CONNECT THEM CONVENIENTLY. 

3) APPLICATIONS 

Looking back to the end of Section 2A, it remains to graft a variational procedure adapted 
to the results of Section 2B. Similarly, as a major result first discovered in specific cases by the 
brothers Cosserat in 1909 and by Weyl in 1916, we shall prove and apply the following key result: 

THE PROCEDURE ONLY DEPENDS ON THE DUAL OF THE SPENCER OPERATOR. 

In order to prove this result, if /g+i, g^+i,/ig+i € H^+i can be composed in such a way that 

g'+l = ffq+l O fq+l = fq+l O hq+l, We get: 

Dg'q+l = fq+l O gfh O jl (gq) O ji {fq) - idq+l = O Djq+l O ji (/,) + DJg+l 

= O /g+1 O Jl (/g) O Jl (^qr) ~ = hq_^_i o Dfq+l O jii^hq) + Dhq+l 

Using the local exactness of the first nonlinear Spencer sequence or ([23], p 219), we may state: 

LEMMA 3.1: For any section fq+i G TZq+i, the finite gauge transformation: 

XqGT*®Rq^X'q = /g+'l « Xg « Jl(/g) + ^/g+l &T*®Rq 
exchanges the solutions of the field equations D'xq = 0. 

LEMMA 3.2: Passing to the limit over the source with hq+i = idq+i + t^g+i + ... for t —>• 0, we 
get an infinitesimal gauge transformation leading to the infinitesimal variation: 

5Xq = D(,q+l + L(jl(Cg+l))Xg 
which does not depend on the parametrization of Xq- 

LEMMA 3.3: Passing to the limit over the target with Xq = Dfq+i and gq+i = idq+i +tr]q+i +..., 
we get the other infinitesimal variation: 

SXq = fq+l ° Dr]q+l O ji{fq) 

which depends on the parametrization of Xq- 
We obtain in particular: 

^x!) = m'^-e) + {edrx':.+xfd.e-x:ie) 

SX;,^ = {^.er-C^) + {f^^rXU+Xlr^^f^+X:iers) 

a result showing the importance of the Spencer operator. In the case of the Killing system Ri with 
g 2 = 0, these variations are exactly the ones provided by the brothers Cosserat in ([7], (49)+ (50), 
p 124, with a printing mistake corrected on p 128), replacing a 3 x 3 skewsymetric matrix by the 
corresponding vector in 
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These two explicit general formulas of the lemma cannot be found somewhere else (The reader 
may compare them to the ones obtained in [14] by means of the so-called ” diagonal ” method that 
cannot be applied to the study of explicit examples). The following unusual difficult proposition 
generalizes well known variational techniques used in continuum mechanics and will be crucially 
used for applications: 


PROPOSITION 3.4: The same variation is obtained whenever rjg = fq+i{^q + XgiO) with 
Xq = Dfq^i, a transformation which only depends on ji{fq) and is invertible if and only if 
det{A) ^ 0. 


Proof: First of all, setting -b Xg(^), we get f for g = 0, a transformation which 

is invertible if and only if det{A) ^ 0. In the nonlinear framework, we have to keep in mind that 
there is no need to vary the object w which is given but only the need to vary the section /g+i as 
we already saw, using ijq G Rq over the target or ^q G Rq over the source. With rjq = /g+i(^g), we 
obtain for example: 

Sfk ^ ^ fk^r 

Sft = + 

J) fk _ £u £v _|_ £u _ £kcr i fk \ £k cr i £k cr 

^ J ij luvJi Jj ' luJij J r >ij ' J ri^j ' J rj^i ' Jrij^ 

and so on. Introducing the formal derivatives di for i = 1, ...,n, we have: 


sf, = cUuvg) = = vtf; + - = frC^ + - + ft.+iX 

We shall denote by tt(? 7 g) = C^iyq,Vq)-^ G ^(^ 9 ) with the corresponding vertical vector 

held, namely: 


avg) = 0 


dx^ 




dy'^ 


ivtiy)yl 


)A 

dyt 


ivuviy)y^y'i + Vuiy)yii) 


dy% 


-b 


However, the standard prolongation of an inhnitesimal change of source coordinates described by 
the horizontal vector held obtained by replacing all the derivatives of ^ by a section ^q G Rq 
over ^ G T, is the vector held: 


d 


d 


d 


Ke,) = + 0^ - (2/fer(a^))TrT " {VriUx) + yUlix) + (x)) 


dx 


dy‘ 


dyt 


dytj 


-b ... 


It can be proved that [b(^)g, b(^'] = b([^g, ^']), V^g, G Rq over the source, with a similar property 
for tt(.) over the target ([23]). However, b(^g) is not a vertical vector field and cannot therefore be 
compared to jl(? 7 g).The solution of this problem explains a strange comment made by Weyl in ([36], 
p 289 -b (78), p 290) and which became a founding stone of classical gauge theory. Indeed, ft 
not a scalar because ft is not a 2-tensor. However, when H = 0, then —Xg is a i?g-connection and 
ft=ft + XtiC is a true scalar that may be set equal to zero in order to obtain ff = —xtiC^ a 
fact explaining why the EM-potential is considered as a connection in quantum mechanics instead 
of using the second order jets ff^ of the conformal system, with a shift by one step in the physical 
interpretation of the Speneer sequence (See [22] for more historical details). 

The main idea is to consider the vertical vector held T(/g)(^) —b(^g) G V{TZq) whenever yq = fq(x). 
Passing to the limit t —>■ 0 in the formula gq o fq = fq o hq, we hrst get g o f = f o h 
f{x) + tri{f{x)) -b ... = f{x + tf{x) + ...). Using the chain rule for derivatives and substituting jets, 
we get successively: 


Sf^ix) = vHfix)) = Cdrf\ 6 ft = Cdrft + ft ft, 6 ft = Cdrft + ftj^l + ftCj + ft^l 
and so on, replacing f^ft^i^ by f^’drft in pq = /g+i(^g) in order to obtain: 


6 ft = vtft + - = emt - ft+u) + ft+iX + - + ftc, 

where the right member only depends on ji(/g) when \ y \ = q. 

Finally, we may write the symbolic formula /g+i(xg) = ji(/g) — fq+i = 7?/g+i G T* 0 V{TZq) in 
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the explicit form: 


+ - + f^,+iXi = Xt. - fUu 

Substituting in the previous formula provides = fq+i{S,q + XgC'?)) we just need to replace q 
by g + 1 in order to achieve the proof. Replacing in the previous variations and using all the for¬ 
mulas involving the Spencer operator and the algebraic bracket that have been already exhibited, 
we let the reader prove as an exercise that we have equivalently: 

X = Xq+i - {Xg+i(-)>?9-ti} 

We obtain in particular: 

X = - xX'') 

) - Xjr + X:i^^r - X';,,i - X^i) 

Checking directly the proposition is not evident even when g = 0 as we have: 

- XdiF = - X - ixx - xix)] 

but cannot be done by hand when g > 1. 

Q.E.D. 

We recall that the linear Spencer sequence for a Lie group of transformations GxX X, which 
essentially depends on the action because infinitesimal generators are needed, is locally isomorphic 
to the linear gauge sequence which does not depend on the action any longer as it is the tensor 
product of the Poincare sequence by the Lie algebra G of G. Accordingly, the main idea will be 
to introduce and compare the three following Lie groups of transformations but other subgroups 
of the conformal group may be considered, like the optical subgroup which is a maximal subgroup 
with 10 parameters, contrary to the Poincare subgroup which is not maximal: 

• The Poincare group of transformations leading to the Killing system R 2 '. 

Glij = “f ^ dr^ijiK) ~ ^ 

^ (L(6)7)y' ^ 4 + lUXl + lUX] - XXr + CdXX) = 0 

• The Weyl group of transformations leading to the system R 2 : 

T ^ dftOij — A(^x')C0ij (^x') 

= (L( 6 ) 7)5 = + xtiXj - XXr + CdXX = 0 

• The conformal group of transformations leading to the conformal Killing system R 2 : 

(-^(Cl)^)u — T “t“ dftOiji^X^ — A(^x')C0ij i^x') 

= 5Xj {^) + {x)X {x)Ar (x) 

where one has to eliminate the arbitrary function A{x) and 1-form Ai{x)dx^ for finding sections, 
replacing the ordinary Lie derivative C{^) by the formal Lie derivative L{£^g), that is replacing 
jqiO by fq when needed. According to the structure of the above Medolaghi equations, it is im¬ 
portant to notice that LI = L{fi)uj G S 2 T* and that P = ^(^ 2)7 G S 2 T* ®T. Moreover, as another 
way to consider the Christoffel symbols, (d,—j) = is a i?i-connection and thus also a 

Ri-connection because Ri C Ri = Ri- 

• We make a few comments on the relationship existing between these systems. 

First of all, when w = (ujij(x) = ujji(x)) is a non-degenerate metric, the corresponding Christof¬ 
fel symbols are 7 = Xj{x) = {x){diUJrj{x) + djUJri{x) — drOJijix)) = 'yji(x)). We have the 
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relations i?i C -Ri = Ri and obtain therefore R 2 = R 2 C pi(Ri), R 2 = pi(Ri), a result 

leading to the strict inclusions R 2 C R 2 C R 2 with respective fiber dimensions 10 < 11 < 15 when 
n = 4 and w is the Minkowski metric with signature (1,1,1, —1). 

Secondly, if we want to deal with geometric objects in both cases, we have to introduce the sym¬ 
metric tensor density ujij = (jJij/\ det{uj) and the second order object 7 ^^- = 7 ^^- — + 

djJri ~ such that I det{uj) |= 1 , 7 ^* = 0 , in such a way that R 2 = {/2 G 112 | = 

a), / 2 ~^( 7 ) = 7 }- It follows that gi is defined by the equations corj^j + = 0 while g 2 

is defined by the equations + dj^ri — = 0 which only depend on u and no 

longer on Cj. Only the first of the three following technical lemmas is known ([21], p 624-628): 

LEMMA 3.5: gi is finite type with 53 = 0 when n > 3. 

Proof: The symbol 53 is defined by the equations — ^{d^SJjrt + Clrt ~ ^ij^^^^rst) = 0- Sum¬ 
ming on k and t, we get = 0. Multiplying by o;*-^ and summing on i and 

j, we get +oj''^Crij = is to say = 0 whenever n>2. Substituting, 

we obtain (n — 2)^G^- = 0 and thus = 0 when n > 3, a result finally leading to ^ = 0 and 
thus g 3 = 0,Vn > 3. In this case, it is important to notice that the third order jets only vanish 
when 7 = 0 locally or, equivalently, when to is locally constant, for example when n = 4 and to is 
the Minkowski metric of space-time. 

Q.E.D. 


LEMMA 3 . 6 : 52 is 2-acyclic when n > 4. 


Proof: As gs+r = 0,Vr > 0, we have only to prove the injectivity of the map 5 in the sequence: 


0 ^ A^T* (g) 52 4 A^T* (g) gi 
and thus to solve the linear system: 


cfc I ek _i_ = n 

Substituting, we get the alternate sum over the cycle, where d is again the Kronecker symbol: 


C{ap^){dU:c.,M + = 0 

Summing on k and i, we get: 




that is to say: 


C{aM){d'lCi,M - = 0 


Summing now on k and a, we get: 


^ ^'i7Crs,a/3) ^ 


Multiplying by and summing on i, we get: 

(n - - co^\dlCs,^, + dl^C.,t 0 ) = 0 

Summing on j and /3, we finally obtain: 

2{n-2)cJ^^C^,J-,=0 ^ Ci, 0 , = O => 4/3^ = 0, Vn>4 

Accordingly, the linear system has the only zero solution and 52 is thus 2-acyclic Vn > 4, a quite 
deep reason for which space-time has formal properties that are not satisfied by space alone. 

Q.E.D. 
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LEMMA 3.7: 52 is 3-acyclic when n > 5. 


Proof: As gs+r = 0,Vr > 0, we have only to prove the injectivity of the map <5 in the sequence: 


0 ^ A^T* A^T* ®T* ®T 

and thus to solve the linear system: 


ffc _ A- = D 

^i/3,'ySa. ' Sz7,(5a/? S.2(5,a^7 

Substituting, we get the alternate sum over the cycle (care to the Kronecker symbol (5): 


Contracting in k and i the previous formula, we get: 


C{ahmc.,M5 = 0 ^ C{ahmtCr^,MS " = 0 

Contracting now in k and a, we get: 


'^^ri,0'yS Crz,7(5/3 ^ri,6'yl3 ^ri^SjS-y ^ (^'i/3Crs,7(5i ^*7^7-8,(5i/3 ^ 


and thus: 


(?r ^')^ri,0'yS ^ (^*/SCr’8,7(5t ^*7^rs,(5/3t rs, ^7t) 0 

that we may transform into: 

(n - + ^Cs,Sl3t + ^S^rs,l3-ft) = ^ 

Contracting in I and 13, we finally obtain: 

2(n-3K^-e,,,-^, = 0 ^ (n-4)e,,^^, = 0 

and g 2 is thus 2-acyclic for n > 5. 


Q.E.D. 


It follows from these lemmas that we are exactly in the same situation as the one met in the 
previous example, with a shift by one in the order of the operators involved. We may thus choose 
Cr = A''T* 'S) R 3 — A''T* (g) R 2 in the Spencer sequence: 

0 —^ e ^ 15 60 ^ 90 60 15 —^ 0 

Each operator Dr is thus induced by the Spencer operator D : 0 i ?3 —>■ A’'T* 0 R 2 and is 

therefore a first order operator with constant coefficients, both with its formal adjoint. Eor later 
computations, the first Spencer operator in the sequence Js^E) T* 0 J 2 {E) A^T* 0 Ji{E) 
can be described by the following images: 




f) 

Si ~ '-'isj 


X 




dif^j - = X 


k 


while the second Spencer operator leads to the identities: 


- a,= 0 , da^j - d,xl, + - xl^ = 0 


Finally, if 22 : E —> E is a linear differential operator of order q, its formal adjoint ad{D) : 
A"T* 0 F* —^ A"T* 0 E* is again a linear differential operator of the same order q that can be 
constructed by contraction with a test row n-form and integration by parts as usual by means of 
the Stokes formula. According to well known properties of the adjoint procedure, if 22i generates 
the CC of 22, then we have ad{V) o ad{Vi) = ad{Vi o 22) = ad(0) = 0 and thus ad{V) is surely 
among the CC of ad{Di) but may not generate them in general. By duality, this remark is at 
the origin of the difficult concept of extension modules in homological algebra and its application 
to the theory of differential modules. It can be proved that such a property does not depend on 
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the differential sequence used, that is one can study 'D^-i and 'Dr in the Janet sequence Vr > 1 
with D — Do or, equivalently Dr and Dr+i in the Spencer sequence, a first highly nontrivial result 
([25], [32]). In the case of the previous systems, as the Poincare sequence is self-adjoint up to 
sign because ad{grad) = —div when X = it follows that ad{Dr-i) generates the CC of ad{Dr) 
while ad{Dr) generates the CC of ad{Dr+i), a second highly nontrivial result (See examples in [27]). 

• We now make a few comments on the relationship existing between these groups. 

As a Lie pseudogroup, the Poincare group is defined by the system TZi C Hi with the n(n-|-l)/2 
equations ujki{y)yiyj = After linearization, (df,— 7 *^) is the only existing symmetric i?i- 

connection for the Killing system i?i C Ji (T) but 7 may also be considered as a geometric object 
of order 2 with well known transition laws. As 32 = 0, : i ?2 Ri is an epimorphism but i?i is 

not involutive and R 2 is involutive whenever the non-degenerate metric ui has constant riemannian 
curvature ([10], [20]). In actual practice, n = 4 and w is the Minkowski metric in the local coor¬ 
dinates (x^, , x'^ = ct). The fact that the Poincare group could have something to do with 

the Galilee group through a kind of limiting deformation procedure with 1/c —>■ 0 is not correct 
because of a few general results on the normalizer T = A^(r) of P in aut{X) which are not so well 
known as their study involves a quite delicate use of the Spencer d-cohomology that we explain 
now (See [29] for more details). 

In 1953 the physicists E. Inonii and E.P. Wigner (1963 Nobel prize) introduced the concept of 
contraction of a Lie algebra by considering the composition law (u, v) — >■ (m -I- u)/(l -I- (uv/c^)) for 
speeds in special relativity (Poincare group) when c is the speed of light, claiming that the limit 
c—^-ooorl/c —>-0 should produce the composition law (u,v) ^ u + v used in classical mechanics 
(Galilee group) ([11])- However, this result is not correct indeed as 1/c ^ 0 has no meaning in¬ 
dependently of the choice of length and time units. Hence, one has to consider the dimensionless 
numbers u = ufc,v = v/c in order to get (u,u) —>■ (u-|-u)/(l-l-uu) with no longer any perturbation 
parameter involved ([18]). Nevertheless, this idea brought the birth of the theory of deformation of 
algebraic structures, culminating in the use of the Chevalley-Eilenberg cohomology of Lie algebras 
([ 6 ], [29]) and one of the first applications of computer algebra in the seventies because a few 
counterexamples can only be found for Lie algebras of dimension >11 and have thus more than 
500 structure constants. Finally, it must also be noticed that the main idea of general relativity 
is to deform the Minkowski metric dx^ + dy^ + dz"^ — c^dt^ of space-time by means of the small 
dimensionless parameter where (j) = GM/r is the gravitational potential at a distance r of a 
central attractive mass M with gravitational constant G. 

It has been the clever discovery of Ernest Vessiot (1865-1952) in 1903 ([35]), still not known 
or even acknowledged today after more than a century (Compare MR0720863 (85m: 12004) to 
MR954613 (90e:58166)), to associate a natural bundle T over X with any Lie pseudogroup P C 
aut{X), both with a section u! of X called geometric object or structure on X as we now explain 
by introducing a copy Y of X and considering the trivial fiber manifold X x Y X. For this 
purpose, Vessiot noticed that any horizontal vector field f = C{x)-^ commutes with any vertical 
vector field rj = on V x AT. Using the chain rule for derivatives up to order q with 

X = x + t^(x) + ... or y = y + tr]{y) -|- ... where t is a small parameter, we may work out the 
respective prolongations at order q on jet coordinates, obtaining therefore the same commutation 
property on H^. As [0, 0] C 0, we may use the Frobenius theorem on the target in order to find 
generating differential invariants {^’'^(yg)} such that = ^'^{yq) whenever y = g{y) € E act¬ 

ing now on the target copy Y of X. Accordingly, prolongations of source transformations exchange 
the differential invariants between themselves, that is any (local) transformation x = ip{x) can 
be lifted to a (local) transformation of the differential invariants between themselves of the form 
u —>■ X{u, jq{(p){x)) allowing to introduce a natural bundle F over X by patching changes of coor¬ 
dinates X = ip{x),u = X{u, jq{ip){x)). A section w of is called a geometric object or structure on 
X and transforms like dj{f{x)) = X{uj{x), jq{f){x)) or simply w = jq{f){uj) whenever y = f(x) is a 
reversible map. This is a way to generalize vectors and tensors {q = 1), connections (g = 2) or even 
higher order objects. As a byproduct we have E = {/ G aut{X)\^i^{jq{f)) = jq{f)~^{uj) = w} as a 
new way to write out the finite Lie equations of E and we may say that E preserves lo. Replacing 
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jq{f) by fq, we also obtain TZq = {fq £ Uq\f~^{uj) = oj}. Coming back to the infinitesimal point 
of view and setting ft = exp{t^) € aut{X),\/f^ G T, we may define the ordinary Lie derivative with 
value in the vector bundle Fq = uj~^{Fo) over X, pull back by w of the vector bundle = y{^) 
over F, by the formula : 

Vf = = C{f)u: = ^ 0 = {^ G T\C{Ou^ = 0} 

We have x ^ x = x + t£^{x) + ... ^ u'^ ^ u'^ = (u) + ... where p = (^i, ...,/x„) is a 

multi-index and we may write down the system of infinitesimal Lie equations in the Medolaghi form: 

L!" = (£(0^)" = -Ll^(u;(x))d^^^ + Cdrcv^ix) = 0 

as a way to state the invariance of the section oj of T. Finally, replacing jq{f) by a section 
fq G Jq(T) over f £ T, we may define Rq C Jq(T) on sections by the purely linear equations: 

= iL{fq)coy = -Lyiu;ix))e^ + Cdrcv^x) = 0 

By analogy with ’’special” and ’’general” relativity, we shall call the given section special and any 
other arbitrary section general. The problem is now to study the formal properties of the linear 
system just obtained with coefficients only depending on ji(w). In particular, if any expression in¬ 
volving uj and its derivatives is a scalar object, it must reduce to a constant whenever F is assumed 
to be transitive and thus cannot be defined by any zero order equation. 

EXAMPLE 3 . 8 : Coming back to the affine and projective examples already presented, we show 
that the Vessiot structure equations may even exist when n = 1. For this, we notice that the only 
generating differential invariant $ = yxx/Vx of the affine case transforms like u = udxf+{dxxf /dxf) 
while the only generating differential invariant = {yxxx/Vx) — ^{VxxIVxy of the projective case 
transforms like v = v{dxfY -f (dxxxf/dxf) - \{dxxf/dxfy when x = f{x). If now 7 is the geo¬ 
metric object of the affine group y = ax+ b and 0 y a = a(x)dx £ T* is a 1-form, we consider the 
object w = (a, 7 ) and get at once one first order and one second order general Medolaghi equations: 

C{f)a = adxf + fdxa = 0, £(C )7 = dxxf + idxf + fdx'^ = 0 

Differentiating the first equation and substituting the second, we get the zero order equation: 

f{adxxCt-2{dxaf F a-fdxa-a^dxl) = Q fdx{^^--) = ti 

a 

and the Vessiot structure equation dxCX — 70 ; = ca^ where c is an arbitrary constant. With a = 
l ,7 = 0 =J>c = 0 we get the translation subgroup y = x Fb while, with a = l/a ;,7 = 0 =J>c =—1 
we get the dilatation subgroup y = ax. Similarly, if n is the geometric object of the projective 
group and we consider the new geometric object u) = ( 7 , v), we get at once one second order and 
one third order general Medolaghi equations: 

^{fh = 9xxf F jdxf F fdxi = 0, C{f)v = dxxxf F 2vdxf F fdxV = 0 

and the only Vessiot structure equation is dx'j — 57 ^ — n = 0, without any structure constant. 

EXAMPLE 3.9: {Riemann structure) If a; = (wy = ujji) G S 2 T* is a metric on a manifold X 
with dim{X) = n such that det{uj) y 0, the Lie pseudogroup of transformations preserving w is 
F = {/ G aut{X)\ji{f)~^(uj) = u)} and is a Lie group with a maximum number of n{n F l)/2 
parameters. A special metric could be the Euclidean metric when n = 1, 2,3 as in elasticity theory 
or the Minkowski metric when n = 4 as in special relativity [18]. The first order general Medolaghi 
equations: 


VLij = {C{f)uj)ij = ojrj{x)di£f F 0Jir{x)dj£f F £fdr0Jij{x) = 0 

are also called classical Killing equations for historical reasons. The main problem is that this 
system is not involutive unless we prolong it to order two by differentiating once the equations. 
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For such a purpose, introducing uj ^ as usual, we may define the Christoff el symbols: 

7*y(a;) = ^UJ^^ix){^^UJrj{x) + djUj„{x) - drUj^ffx)) = j^i{x) 

This is a new geometric object of order 2 providing the Levi-Civita isomorphism ji(w) = (w, duj) ~ 
(u;, 7 ) of affine bundles and allowing to obtain the second order general Medolaghi equations: 

= (£( 07 )?, = - lUx)dre + CdrlU^) = 0 

Surprisingly, the following expression, called Riemann tensor: 

Puj{x) = ^^-flj{x) - 5,7i?(a;) + 

is still a first order geometric object and even a 4-tensor with n^{n^ — 1)/12 independent compo¬ 
nents satisfying the purely algebraic relations : 

P% + P%1 + Pjli = 0: ’^rlPkij + ^krPuj = 0 

Accordingly, the IC must express that the new first order equations = {£{^)p)f^j = 0 are only 
linear combinations of the previous ones and we get the Vessiot strueture equations: 

Phj{x) = c{6^u}ij{x) - 6^u}ii{x)) 

with the only structure constant c describing the constant Riemannian curvature condition of 
Eisenhart ([10], [20], [22], [23]). One can proceed similarly for the conformal Killing system 
C{f)uj = A(x)uj and obtain that the Weyl tensor must vanish, without any structure constant 
([20], p 132). Though this result, first found by the author of this paper as early as in 1978 ([20]) 
is still not acknowledged, there is no conceptual difference at all between the unique structure con¬ 
stant c appearing in this example and the previous one. Moreover, the structure constants have in 
general nothing to do with the structure constants of any Lie algebra. 

More generally, any generating set {4>'^} of differential invariants must satisfy quasi-linear CC 
of the symbolic form v = I{ui) = A{u)ux + B{u) = 0 where ui = {u,Ux), allowing to define an 
affine subfibered manifold Bi C Ji(-R) over IF and a natural bundle = Ji(F)/Bi over F with 
local coordinates {x,u,v). The Vessiot structure equations I{ui) = c(u) are defined by an equiv- 
ariant section c : F ^ Fi : (x,u) ^ (x, u,v = c(u)) depending, as we just saw, on a finite number 
of constants (See [20] and [23] for details and other examples). The form of the Vessiot structure 
equations is invariant under any change of local coordinates. The following result, already known 
to Vessiot in 1903 ([35], p 445), is still ignored today. For this, let us consider two sections w and u) 
of F giving rise, through the corresponding Medolaghi equations, to the systems Rq and Rq. We 
define the equivalence relation: 

DEFINITION 3.10: ui ^ oo ^ Rq = Rq. 

The following result is not evident at all ([20], [29]): 

PROPOSITION 3.11: uj is obtained from w by a Lie group of transformations acting on the 
fibers of F, namely the reciprocal of the Lie group of transformations describing the natural struc¬ 
ture of F. These finite transformations of the form u = g{u, a) will be called label transformations 
and the number of parameters a is < dim{Jq{T)/Rq) = dim(F(f). 

COROLLARY 3.12: Any finite label transformation u = g{u,a) induces a finite transformation 
c = h(c, a) and we say that w ~ w —>■ c ~ c. 

DEFINITION 3.13: The normalizer f = N{T) of F in aut{X) is the biggest Lie pseudogroup 
in which F is normal, that is (roughly) N{T) = f = {/ € aut{X)\f o f o f~^ G r,V/ G F} and we 
write F <1 7V(F) C aut{X). 
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Of course, N{Q) = {rj G ?"|[C)''7] C 0,V^ G 0} will play the part of a Lie algebra for A^(r) 
exactly like 0 did for L. However, N{T) may have many components different from the con¬ 
nected component of the identity. For example, when n = 2 and F is defined by the system 
{j /2 = = 0}, then fV(F) is defined by the system {yl = = 0} U {yl = 0,y22 = 0} as it 

contains the permutation y^ = x^,y'^ = . In actual practice, F = {/ G aut{X)\uj = jq{f)~^{(jj) = 

g{uj,a), h(c,a) = c} is defined by the system TZg+i = {fq+i G ng+i|/q+i(i?g) = Rq} with lineariza¬ 
tion Rq^l — ^^q^l\L(^^q^l^Tjq G Rqj^Tjq G Rq'^ ^ that iS tO Say G Rq 

{^g+i,? 7 q+i} G Rq. Accordingly, the system Rq+i defining 0 = N{<d) can be obtained by purely 
algebraic techniques from the system defining 0. We have ([20], p 390; [21], p 715; [22], p 548;[29]): 

PROPOSITION 3.13: If Rq is formally integrable and gq is 2-acyclic, then Rq+i is formally 
integrable with gq+i = gq+i- 

EXAMPLE 3.15: In the previous Example with oj = ( 0 , 7 ), we obtain by substraction Co ~ 
w 77 (d = aa ,7 = 7 - 1 - ba) c = 7 C — The condition c = c provides b = {1 — a)c, that 
is to say & = 0 if c = 0 and & = a—lifc = —1. Hence, in both cases the corresponding Lie 
pseudogroup is of codimension 1 in its normalizer. Indeed, the normalizer oi y = ax is y = ax^ 
while the normalizer oiy = x + b is y = ax+ b with different meanings for the constants a and 
b. Similarly, in the case of the Riemann structure, we let the reader prove as an exercise that 
w~w77w = aa;—>- 0=^0 because 7 = 7 . Accordingly, the corresponding Lie pseudogroup is of 
codimension zero in its normalizer if c yf 0 and of codimension 1 if c = 0 , a result explaining why 
the normalizer of the Poincare group is the Weyl group, obtained by adding a unique dilatation for 
space and time, contrary to the Galilee group which is of codimension 2 in its normalizer, obtained 
by adding separate dilatations for space and time ([22], [29]). We invite the reader to treat similarly 
the examples provided in the first section in order to understand how tricky are the computations 
involved or to look at the example fully treated in ([21], p726). 


• We now study each group separately, in relation with applications. 


EXAMPLE 3.16: (Poincare group) Changing slightly the notations while restricting for sim¬ 
plicity the formulas to the plane with n = 2 and local coordinates {x^,x'^) instead of space 
with n = 3 and local coordinates {x^,x'^,x^) or space-time with n = 4 and local coordinates 
{x^, x"^, x^, x'^ = ct), we may copy the equations (12) of ([7], p 14) and (12') of ([7], p 19) side by 
side in the following way: 


/ dF 


, dG 

\ 


dH ( dx 


i ax ax IT' 

+ - - 17 ^ 


L 

M 

N 


dF 

ds 


dG 

ds 


[ £iH + x^G - x^F) 


L 

M 

N + x^M — x^L 


We notice that the left members of the equations on the right hand side are only made by the 
derivative of an expression with respect to the curvilinear abscissa s along the curve considered in 
the plane with local coordinates (a;^(s), a;^(s)). Equivalently, we may use the linear transformations 
(F, G, H) -G {F' = F,G' = G,H' = H + x^G - x^F) and (L, M, N) -)■ {V = L,M' = M, N' = 
N + x^M — x'^L) with the same underlying 3x3 matrix of full rank 3, namely: 


/1 

0 —x"^ \ 

/ 1 

0 —x"^ \ 

{F',G',H') = {F,G,H)\ 0 

1 

x^ 

, {L',M',N') = {L,M,N)\ 0 

1 

x^ 

V 0 

0 

1 / 

V 0 

0 

1 / 


but this result is not intrinsic at all and just looks like a pure coineidenee. It is important to notice 
that, while these formulas have been exhibited in the study of the (static) deformation theory of 
a line (Chapter H of [7], p 14 and 19), similar formulas also exist in the study of the (static) 
deformation theory of a surface (Chapter HI of [7], p 76 and 91) and in the study of the (static) 
deformation theory of a medium (Chapter IV of [7], p 137 and 140). We shall not insist on these 
points which have already been treated elsewhere with full details ([ 22 ], [28]) and that we have 
recovered in this paper by means of other methods, but invite the reader to look at the amount 
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of calculations provided by the brothers E. and F. Cosserat. However, in order to establish a link 
between this example and the use of the Spencer operator, we now consider the Killing system 
for n = 2 and the euclidean metric. The dual of the Spencer operator is provided by the integra¬ 
tion by parts of the contraction 2 -form while raising or lowering the indices by means of the metric, : 

— ^14) -I- — C2,l) + Cr^’^( 92^1 — ^1,2) + — ^2,2) + — 0 ) 

provided that ^14 = 0 ,^i ^2 + ^ 2,1 = 0 ,^ 2,2 = 0. Integrating by parts, the factors of ^i,'C 2 ,?i ,2 
furnishes (up to sign) the Cosserat equations where, of course, may be different from 

In arbitrary dimension, one should get similarly ([7], p 137 for n = 3, p 167 for n = 4) (See [28], 
Remark 7.1, p 25 for more details): 

= f, dr - cr^’* = Vi < j 

As a byproduct, we obtain: 

— xV^’'’) = -I- a;^/* — x^f^, VI < i < j < n 

that is exactly the equation used in continuum mechanics in order to study the torsor equilibrium 
bringing the symmetry of the stress tensor when /i = 0 and to = 0 , where the left member is the 
Stokes formula applied to the total surface density of momentum while the right member is the 
total volume density of momentum (See [22], [24], [27] for more details and compare to [34]). 
With the infinitesimal generators 9i = di, 82 = 92 , 9 ^ = x^d 2 — x^di, setting = X^d^j^9r with 
^^{x) & G, we have: 

= X^+ x^X^,^l = 

and we find back exactly the 3x3 matrix with full rank already exhibited. 


In fact, our purpose is quite different now though it is also based on the combined use of group 
theory and the Spencer operator. The idea is to notice that the brothers are always dealing with 
the same group of rigid motions because the lines, surfaces or media they consider are all supposed 
to be in the same 3-dimensional background/surrounding space which is acted on by the group of 
rigid motions, namely a group with 6 parameters (3 translations + 3 rotations). In 1909 it should 
have been strictly impossible for the two brothers to extend their approach to bigger groups, in 
particular to include the only additional dilatation of the Weyl group that will provide the virial 
theorem and, a fortiori, the elations of the conformal group considered later on by H.Weyl. In 
order to emphasize the reason for using Lie equations, we now provide the explicit form of the n 
Hnite elations and their infinitesimal counterpart, namely: 


y = 1 _ 2 {bx) + b'^x'^ => ^s = -^x'^S^dr + ujstx*x^dr dr9l = nuistx*, yi<r,s,t<n 
where the underlying metric is used for the scalar products a;^, bx, involved. 


EXAMPLE 3.17:{Weyl group) We may rewrite the infinitesimal Lie equations in the local form: 

0 < Z, j, fc < n 

The contraction form should be complemented by the terms 12 ^ (dr^} — 0) and the integration by 
parts provides the additional dual equation dri 2 ^ + = v because = ... = ^” = ^^r (Compare 

to (74) in [36], p 288). As a byproduct, we get: 

dr(i 2 ’’ + x"a[) = V + x^'fi 

that is exactly the virial equation already presented for the symmetric stress used in continuum 
mechanics and gas dynamics, where the left member is the Stokes formula applied to the total 
surface density of virial while the right member is the total volume density of virial. 
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Introducing the additional infinitesimal generator 04 = x^di, we now get: 


and obtain the 4x4 matrix of rank 4: 


C2 _ 

A2 


+ x'^) 

/ 1 

0 

—x'^ 

x^ \ 

0 

1 

x^ 

x^ 

0 

0 

1 

0 

V 0 

0 

0 

1 / 


describing the linear map X x Q ^ Ri with dim{Ri) = 4 when X = . 


EXAMPLE 3.18: {Conformal group) First of all, we explain the confusion done by Weyl in ([36]) 
between natural bundles and jet bundles, recalling that both bundles have only been introduced 
fifty years later but that the formula of Weyl that we shall consider has been one of the key ingre¬ 
dients of gauge theory, also fifty years later but for a quite different reason (See ([22]), Chapter 5, 
p 321-343 for historical comments). Indeed, considering uj and 7 as geometric objects, we obtain 
at once the formulas: 


= a(x)u 


Iri — Iri + 




Though looking like the key formula (69)in ([36], p 286), this transformation is quite different 
because the sign is not coherent and the second object has nothing to do with a 1-form. More¬ 
over, if we use n = 2 and set C{^)uj = Auj for the standard euclidean metric, we should have 
(3ii + 822 )A = 0, contrary to the assumption that A is arbitrary which is only agreeing with the 
jet formulas: 


L(a)cc = Acc ^ 2{C + = nA, {L{^ 2 l)U=nA, € 4 

Now, if we make a change of coordinates x = (p{x) on a function a € A'^T*, we get: 


a{(p{x)) = a(x) 


da dipl 


da 

dx^ 


dxi dx^ 

We obtain therefore an isomorphism Ji(A°T*) ~ A^T*XxT*, a result leading to the following 
commutative diagram: 


i ?2 

ID 

T* O Ri 


i ?2 
i D 

T* O Ri 


JiiA^T*) 
i D 


where the rows are exact by counting the dimensions. The operator on the right is D : (^A, Ai) —y 
{^diA— Ai) and is induced by the central Spencer operator, a result that could not have been even 
imagined by Weyl and followers. 

Though striking it may loo like, this result provides a good transition towards the conformal ori¬ 
gin of electromagnetism. The nonlinear aspect has been already presented in ([22], [23], [31]) 
and we restrict our study to the linear framework. A first problem to solve is to construct vec¬ 
tor bundles from the various components of the image of Di. For this purpose, let us introduce 
(Bti = ^4 + 4.41) G r* 0 T* 0 T with {Bf, = Bf) G T* and (Sf^-, = 4^ + 7 ^- 
lijX^^i+ X'^^dr^i^) G T*®S 2 T*®T with {BA . — BA ,. = 4 ) g A^T*, We obtain from the relations 
dilfj = djjA and the previous identities: 


p. . _ Dr _ or 

nj ^ri,j ^rj,i 


-yr 


_ yr 

- d,x: 


■Yr 


4f,4 - YrsXY + Xf^drYs^ ” Xf.drYsj 


da;., - djx;, + j;axx _ + 74 , 71 . - x;d,j, 

^^(X;, + YrsX;,) - d,iX;, + YrsX^ ^ 

diBj - djBi 


s,i) 
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Now, we have: 


B, = {diCr-Cr,)+Yrs{dii’^ -^t) 

= diC + + C®5s7ri - it-A 

= d.iC + YrsO-nA, 

= n{\diA-Ai) 

and we finally get Fij = n{djAi — diAj), a result fully solving the dream of Weyl. Of course, when 
n = 4 and to is the Minkowski metric, then we have 7 = 0 in actual practice and the previous 
formulas become particularly simple. 

As Cr = A^T* (g) i?2 C A^T* ®R 2 = Cr and R 2 /R 2 ^ T*, we get Cr/Cr ^ A^T* 0 T* and the 
conformal Spencer sequence projects onto the sequence T* —>• T* (g)T* —>• A^T*(g)T* —>• .... Finally, 

the Spencer sequence projects with a shift by one step onto the Poincare sequence T* —>■ A^T* —>■ 
A^T* —>■... obtained by applying the Spencer map 6 , because these two sequences are only made 
by first order involutive operators and the successive projections can therefore be constructed in¬ 
ductively. The short exact sequence 0 —>■ S 2 T* -A T* (g) T* -A A^T* —0 has already been used in 
([22], [23], [28], [30]) for exhibiting the Ricci tensor and the above result brings for the first time 
a conformal fink between electromagnetism and gravitation by using second order jets. 

As for duality, using standard notations, we have the possible additional terms: 

... + - C„) + E . - djfA) 

Developping the sum, we get: 

... - J^Cl - - - J^Cn + - d 2 Cl) + ... + ” dnCn) + - 

and the integration by part provides therefore the equations: 

djr^ -T = w^ 

The link with the virial theorem is provided by the formulas fA = = ... = nfA and 

F = nj"^. Accordingly, when the second members of the inductions equations vanish, we have 
diJ^ = dijF^^ = 0 crjl = 0 in a coherent way with a well known property of the so-called 
impulsion-energy tensor in electromagnetism. 

We have therefore obtained the following crucial theorem and striking corollary: 

THEOREM 3.19: In the Spencer sequence for the conformal Killing system with n = 4, the 
field equations Ci ^ C 2 projects onto the first set of Maxwell equations A^T* 
the correponding potential parametrization Cq —b- Ci projects onto the usual parametrization 
A^T* —^ A^T* by the electromagnetic 4-potential, according to the following commutative dia¬ 
gram: 


A°T* (g) R 2 


0 


^1. 


A^T* (g) R 2 

i 

A^T* 

i 

0 


Dn 


A^T* (g) i?2 

i 

i 

0 


By duality, the second set of Maxwell equations A” A” ^T* is induced by the induction 

equations a”T* (g) Ci ^ (j* g similar property holds for the corresponding pseudo¬ 

potential parametrizations, according to the following commutative diagram: 


0 


0 


0 

i 


i 


i 

y\^n— 

ad{d) 

y\^n—2^* 

ad{d) 


i 


i 


i 

A"T* (g) Rl 

ad{D \) 

A^-iy* (g) R* 

arf(D2) 

^n-2T* 0 fll 
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Similar comments can be done for the Clausius and Cosserat equations because i?2 C -R2 C i?2- 


REMARK 3.20: The key formulas (76) in ([36], p 289) are based on a confusion between the Janet 
and Spencer sequences. Indeed, using only components of the Spencer operator when ^2 G R27 we 
have on one side in S 2 T*: 


UJrjX^^ + UJirXj 


-iD+OJiridjC -Cj) 

(uJrjd^C + (^irdjC + Cdri^rj) “ f (C + 


Similarly, we have on the other side in T*: 


x: 


.Xf 


a*(e + 7?.^) - (^i + lUt + 

n{\diA - Ai) 


In the nonlinear framework, using the variational formulas already establihed at the beginning of 
this section, it follows that: 


(tti = Xr.i + TrsXfi) G r* ^ 5a^ = n{-d^A - Ai) + ^ra^ + - nx'^.As) 

a formula that cannot be fully identified with (76) because A^ ^ 0. 


REMARK 3.21: As another confusion, we revisit a basic result of classical gauge theory. First 
of all, we recall that the classical lagrangian of a free particle of mass m and charge e in an EM 
field F = dA is L{t,x,x) = ^mcOijFx^ + eFAi and let the reader check that the corresponding 
Euler-Lagrange equations are ^ = v A B where v = (i*) and the right member is the Lorentz force. 
Introducing the momentum pi = muJirxA + eAi and substituting, we obtain easily the hamiltonian 
F[{t,x,p) = 2^w*'^ (pi — eAi){pj — eAj) which is obtained from the hamiltonian H = ■^oj^^piPj of 
the uncharged particle by the transformation pi —> pi — eAi. A main idea of gauge theory has 
been to transform pi to —ihdi according to the correspondence principle of quantum mechanics 
and Ai to a connection. Quite contrary to this point of view, we have: 


COROLLARY 3.22: The transformation pi —>■ pi — Ai is provided by the dual action of the 
second order jets of the conformal system and has a purely group-theoretical origin. 


Proof: The nonlinear system 7^2 C 112 of conformal finite Lie equations is: 


Wfci/f/j = a{x)u}^j ^ g^{fij F lisflfj) = 1% + ^faj{x) -f S’)ai{x) - uj^ju'^^arix) 

Using the formulas of the last proposition describing the transformation 771 = /2(Ci) with /2 G 7^2 
and ^1,771 G Ri, we obtain by contraction 77^ = fiC^Vk = C + noi^^. The inverse transformation 
allowing to describe R* is = Vk~ if we set a* = fibk as in ([23], p 448). Hence, 

according to the variational formula 5xi = o 77772 o ji(/i) of the first lemma of this section, if 
the dual finite field of ix%Xr,i) <^xer the source is (A^*, A*), we obtain for the dual field {yi^ 
over the target the dual formulas replacing the ones used for the adjoint operator: 


A/a,/'= - nbiX^d,f\ Ay’^ = 

This result can also be found by a direct computation showing that: 




The action of the second order jets only is yf —^ yf — naiy^ if we set f{x) = x, fb 
corollary follows when fc = 4 because y^ = 77^ is the charge density. 

Q.E.D. 


Sf. The 


REMARK 3.23: As we have already seen, the Helmholtz analogy establishes a parallel between 
analytical mechanics and thermodynamics by setting T = q and L = —F, that is by trying to 
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describe the absolute temperature as a classical ” field ” that should be obtained from the deriva¬ 
tive of a ”potential ”, in order to bring the possibility of an integration by part. Accordingly, it 
should not be just a scalar quantity but should transform like a time derivative and we recognize 
the definition of a jet. On the other side, it follows from the preceding results that the Poincare 
group is of codimension 1 in its normalizer which is the Weyl group, the ” difference ” being the 
unique space-time dilatation that must be added and which is not accessible to intuition as it 
must be at the same time a space dilatation and a time dilatation, that is a concept where space 
cannot be separated from time, contrary to the standard examples of thermostatics that we have 
provided. Also, we must not forget that the measure of V has to do with the translation of a 
piston ring while the measure of T has similarly to do with the translation of a liquid or a solid in 
a thermometer. The link with thermodynamics will be obtained by defining as usual the so-called 

” dimensionless speed ” v^/c = 84 , and the ” normalized speed ” u* = (u*/c)/-y /1 — ^ over the 
target with ujki{y)u'^u^ = — 1 u^-^{u}rkiy)u^) = 0. Also, introducing = iOrs{x)A\A‘t while 

setting 0 ^(/(a;)) = l/a(a;) where 6 is a gauging of the dilatation subgroup over the target as in 
section 2 , we may introduce: 


I 0 = ff^A^ = e^i - ^ 84 /^ = 00 

I p = ff-A44ldet[A) = (l/e^jy/f - T^ 8 {x^,x^,x^)l 8 [y^,y^,y'^) = p/O^ 

It follows that pvfi = 3194 /^ and we obtain therefore the identity 8 {pu^)/dy^ = 0 over the target 
from the more general identity 8 {^ 8 if^)/ 8 y^ = 0 over the target ([25], Lemma 3.94, p 490). We 
have thus a reason for introducing the constraints A\ = A'^ = af = 0 on the moving frame like in 

([7], §28, §45), obtaining the so-called temperature vector /!(x) = by gauging the Lorentz 

subgroup. For a fluid at rest, we have 0 — 084 /'^ and we cannot separate 0 from 84 /"^ in the non¬ 
linear description of the Spencer operator because xo = A — id, explaining therefore the Helmholtz 
analogy because the evaluation y^ = on time does not commute with the integration by part 
on time. Finally, introducing a density of action r<;(A) over the source, we get for the variation: 


SW = 6 J w{A)dx = J ^SAldx = J - rif)8^fdx = J yli^ - Vi)dy 

where yl = \ 9 k^^ 8 iff over the target. With a density w{p, 0) = w{p, 0, 0), we get: 

1 1 o„.. 

yl = —7r(5[ -I- LUrkU^U^) — aOJrkU^u’' 


1 ^8w 


1 ~ 8 w 
a = —0 —^ 

A 80 


For example, we may set —w = det{A) = ^a{0/p) = ^q; 0^A for the black body and recover 
the pressure tt = |q;0 ^. 

Finally, using a i?i-connection (d, —7) with 7 = 0 when n = 4, we may introduce = (^ = 
(0,0, 0,1), fff = 0) G i?i C .Ri C Ji{T) over the source and get y^ifix)) = 94/''(x), p^(/(x))/“(x) = 
84 /^(x) over the target. The following tricky proposition revisits the first principle of relativistic 
thermodynamics in this new framework while taking into account the Helmholtz analogy: 


PROPOSITION 3.24: We have the formula: 

8 y! 


kk 


8 ,~ 8 w 


8 y 


,y^ + ylyf = pu^^^{6^^-w) 


Proof: First of all, as = 0u^, 0 = pA and u^{dl + (jJrkU^u^) = 0, we have: 

+ UJrkVd'ff)) = -Tr{Sl + UJrku’^fffiu'^^ + 0^) 


= -TT{5l+UJrkU'"u^)0 — 


k 


= 

dy’- 

= 

= wpp'w 


dy^ 
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Similarly, we have: 




= —pu 




fn aw \ I dw _.,i 

+ IWPU 




Finally, multiplying the conformal Killing equations for rji over the source by u\ we obtain 
UrkU^u^Tli = —^Vk — we have thus for n = 4: 


^ dw „ _ dw ~ dw . , „ X 

+ 0 —^ = A(—Stt + a) 
(90 dp 80 


k \ k dw 1 dQ 

ykVi = -(-Stt + = --^pu 


The proposition follows from the fact that dw = ^dp+ ^dO + fgrf© along the trajectory. The 
extension to continuum mechanics could be done by using the 6 quantities A*^ = Aij — for 

i,j = 1,2,3 which do not contain Af for i = 1,2, 3,4, in place of p and we have det{A*-) = 1/p^. 
(Compare to [17]). 


Q.E.D. 

From the infinitesimal point of view, we may use (up to sign) (( 9 i^^+ 92 ^^+ 935 ^) —(^J+^l+^f) in 
place of p, did~^‘i place of 9 and — in place of Xr,i- Identifying the speed with a Lorentz 
gauging, we may introduce the constraint 94 ^* — ^4 = 0 in order to have 944 ^* — ^44 = 94^4 — ^44 
and we may consider by substruction the components 9i^| — 944 ^* because = Ai for 

i = 1,2,3, a result explaining why it is possible to combine the gradient of temperatue with the 
(small) dimensionless object obtained when dividing the acceleration by the square of the speed 
of light as in ([9], (37), p 922). Also, considering the coupling bilinear term ( 94 ^* — QC){did — CU 
as in ([24]), we understand how it is possible to introduce a spatial heat flow in the components 
a\ of the stress ct for i = 1,2,3, contrary to the phenomenological approach of ([9], (24)+(37)). 

EXAMPLE 3.25: {Projective group) With n = 1, let us consider the projective transformations 
of the real line. For this, we may introduce the third order system R 3 C JsiT) of infinitesimal Lie 
equations defined by ^xxx = 0 and we let the reader check easily that [i? 3 , i? 3 ] C R 3 . We may then 
exhibit a basis made by the three infinitesimal generators { 0 i = dx^ 0 i = xdx^Os = \x^dx} while 
introducing = X^d^Ox leading to ^ = X^Ot, ^x = X^dxOr, ^xx = X^dxxOr along with the following 
corresponding linear transformation: 

U \ / 1 X 1x2 \ / Ai \ 

(L, M, A) e, = (L, M, A) 0 1 X 

\^xx ) V 0 0 i j \ X^ j 

bringing the relations: 

X = dx^-^x = {dxX")9r, Y = dx^x - ixx = {dxX^)dx 0 T, Z = dxixx - ixxx = {dxy )dxx 0 T 
Finally, it just remains to integrate by parts the expression/contraction: 


FX + GY + HZ = F{dx^ — ^x) + G{dx^x ~ ^xx) + H{dx^xx ~ ^xxx) 


while taking into account the fact that ^xxx = 0 in order to find the dual Cosserat equations: 


' e 

^ ^x 
. ^xx 


dxF = L 
dxG + F = M 
dxH + G = N 


( dx{F) 

dx{G + xF) 


= L 

= M + xL 


y dx{H + xG + ^x'^F) = N + xM + ^x'^L 


involving the formal adjoint of the first Spencer operator and the above linear transformation act¬ 
ing on both sides of the equations. The study of the conformal group is quite similar to that of the 
projective group because the symbol at order 3 is equal to zero in both cases. We may therefore 
just replace dr^'' - di^ - di^j - 0 by dx^ - ^x, dx^x - ^xx, dx^xx - 0 respectively. 
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4) CONCLUSION 


Considering a Lie group of transformations as a Lie pseudogroup of transformations, we have 
revisited in this new framework the mathematical foundations of both thermodynamics and gauge 
theory. As a byproduct, we have proved that the methods known for Lie groups cannot be adapted 
to Lie pseudogroups and that the two approaches are thus not compatible on the purely math¬ 
ematical level. In particular, the electromagnetic field, which is a 2-form with value in the Lie 
algebra of the unitary group U{1) according to classical gauge theory, becomes part of a 1-form 
with value in a Lie algebroid in the new conformal approach. More generally, shifting hy one step 
the interpretation of the differential sequences involved, the ’’field” is no longer a 2-form with value 
in a Lie algebra but must be a 1-form with value in a Lie algebroid. Meanwhile, we have proved 
that the use of Lie equations allows to avoid any explicit description of the action of the under¬ 
lying group, a fact particularly useful for the elations of the conformal group. However, a main 
problem is that the formal methods developped by Spencer and coworkers around 1970 are still 
not acknowledged by physicists and we don’t even speak about the Vessiot structure equations for 
pseudogroups, not even acknowledged by mathematicians after more than a century. Finally, as a 
very striking fact with deep roots in homological algebra, the Clausius/Cosserat /Maxwell/Weyl 
equations can be parametrized, contrary to Einstein equations. We hope this paper will open new 
trends for future theoretical physics, based on the use of new differential geometric methods. 
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